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FOREWORD 


This  Project,  (6906)  "Nuclear  Weapon  Effects  on  Space  Vehicle*;1  and  Task,  (690601)  "Deter* 
mination  of  High  Altitude  Nuclear  Weapon  Effects  on  Space  Vehicles,  "  are  a  part  of  the  Air  Force 
System*  Commend  applied  research  program  710A,"  NUCLEAR  WEAPONS  EFFECTS.  " 

The  study  was  initiated  by  the  Flight  Dynamics  Laboratory,  WWRMD,  under  Project  No.  6906, 
Task  No.  690601.  The  project  officer  was  Mr.  L.  E.  Gilbert.  The  investigation  was  conducted  by 
the  Southwest  Research  Institute  during  the  period  from  1  March  1961  through  31  March  1962. 
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ABSTRACT 

Procedure*  for  analytically  predicting  the  reiponee  of  mieeile  bodie*  to  bleat  loading*  are 
pretented.  The  investigation  involve*  the  behavior  of  cylindrical  shell*  (with  various  end-closure*) 
and  circular,  flat  plate*.  The  numerical  result*  obtained  from  the  analytical  method*  compare 
favorably  with  the  experimental  data  acquired  during  the  study. 
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1.  INTRODUCTION 


In  evaluating  the  vulnerability  of  miesile  eysteme  to  nuclear  antimissile  weapons,  one 
aspect  of  considerable  importance  is  the  response  of  the  vehicle's  structure.  The  analytical  (and 
to  a  lesser  extent,  the  experimental)  information  presented  in  this  report  is  directed  towards 
developing  procedures  for  predicting  the  response  of  missile  bodies  as  represented  by  cylindrical 
shells  with  various  types  of  end- closures. 

Within  the  concept  of  structural  response,  the  most  significant  parameters  which  con¬ 
stitute  a  measurement  of  the  damage  a  missile  will  sustain  as  the  result  of  blast  loadings  are  the 
deformations  and  accelerations.  Consequently,  each  of  the  analytical  procedures  set  out  to  ini¬ 
tially  define  the  time  dependent  variation  of  the  displacement  components.  Once  this  is  accom¬ 
plished,  it  becomes  a  simple  matter  to  obtain  accelerations  as  functions  of  time.  For  the  particular 
shock  sensitive  components  associated  with  each  missile  system,  this  would  provide  an  indication  of 
the  critical  acceleration  forces  to  be  expected. 

It  is  analytically  expeditious  to  subdivide  the  overall  missile  structure  into  two  parts,  each 
representing  different  vulnerability  aspects.  The  first  is  to  consider  the  cylinders'  end-closures 
only  insofar  as  their  restraint  contributes  to  deformational  response  of  the  shell.  Here,  it  is 
implied  that  the  missile  shell  itself  is  most  vulnerable  to  nuclear  weapon  effects.  The  other 
approach  is  to  consider  the  cylinder  end-closures  themselves.  For  example,  a  flat,  circular  plate 
would  require  careful  analysis  if  it  can  be  shown  its  response  is  considerably  more  severe  than  that 
experienced  by  the  shell. 

Accordingly,  in  Section  2,  comparatively  straightforward  analytical  procedures  for 
acquiring  displacements  that  are  reasonably  accurate  and  not  unnecessarily  laborious  are  pre¬ 
sented.  For  the  cylindrical  shell,  both  the  small  and  large  deflection,  linear-elastic  theories  are 
used.  In  the  small  deflection  approach,  the  effect  of  end-closures  (simple,  fixed  and  elastic  sup¬ 
port)  are  considered.  Two  treatments  of  the  circular,  flat  plate  are  presented:  one  entails  the 
dynamic  response  of  circular  plates  at  large  amplitudes  and  the  other  the  plastic  collapse  of  cir¬ 
cular  plates  under  blast  loadings.  Each  analytical  procedure  is  carried  out  in  sufficient  detail  so 
as  to  present  results  which  are  immediately  applicable  to  specific  problems. 

In  Section  3,  extensive  numerical  computations  are  presented  for  the  analytical  procedure 
based  on  the  linear- elastic,  small  displacement  shell  theory.  In  evaluating  the  analytical  pro¬ 
cedures  for  predicting  the  response  of  cylindrical  shell  structures  to  blast  loadings,  there  are 
several  logical,  as  well  as  practical,  reasons  for  initially  investigating  the  adequacy  of  the  simplest 
method.  Foremost  among  these  reasons  is  the  fact  that  these  methods  provide  numerical  solutions 
for  a  wide  variety  of  variables  with  a  minimum  expenditure  in  computational  effort.  The  extent  to 
which  various  load  arjd  time  parameters  and  certain  approximations  influence  the  numerical 
answers  and  cause  these  predicted  answers  to  approach  or  deviate  from  the  experimental  results 
can  be  quickly  determined. 

Moreover,  once  these  numerical  solutions  are  available,  they  also  become  of  value  in  form¬ 
ing  a  foundation  for  subsequent,  more  complex  analytical  procedures.  This  approach  has  been 
adopted  in  the  numerical  analysis  of  the  cylindrical  shells.  In  Section  3  are  the  results  of  the  cal¬ 
culations  using  analytical  methods  based  on  linear-elastic,  small  displacement  theory.  Included 
are  comparisons  of  various  combinations  of  overpressure  and  positive  phase  durations. 

The  information  and  data  obtained  in  the  limited  experimental  portion  of  the  program  are 
presented  in  Appendices  I  and  II.  These  data  (involving  the  radial  displacements)  are  compared 
with  the  analytical  results  in  Section  4.  ”™ 
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2.  ANALYTICAL  PROCEDURES 


2.1  Re»pon»e  of  Cylindrical  Shell*  to  Impultive  Loadings 

2.1.1  Linear-Elastic,  Small- Deflection  Theory 

The  development  of  the  necessary  equation*  for  the  analysis  of  cylindrical  shells  begins 
with  the  displacement  (Equation(I. 12)]  and  frequency  {Equation(I.  13)]  equations  in  Appendix  I.  We 
assume  that  the  cylinder's  end-closures  consist  of  elastic  media  which  offer  unequal,  elastic 
restraint  against  rotation.  For  generality,  we  define  these  elastic  restraints  as  Sj  at  x  =  0  and 
S2  at  x  =  L  where  S  is  the  stiffness  per  unit  length  of  the  medium  or  the  moment  required  to  rotate 
a  unit  length  of  the  medium  through  an  angle  of  one  radian. 

The  boundary  conditions  are  taken  as 
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The  first  two  equations  denote  the  condition  of  continuity  between  the  shell  and  the  end-closures. 

The  condition  that  the  radial  displacements  are  zero  at  the  shell's  boundaries  appears  reasonable, 
particularly  for  end-closures  (such  as  flat  plates)  which  offer  a  high  degree  of  transverse  restraint. 
Since 


Equations  (2. 1)  may  be  rewritten  as  follows 
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where 


R  = 


SL 


Utilizing  Equation!  (2.  2)  in  Equation#  (1. 12)  (Appendix  I)  we  obtain  the  following  eet  of 
linear,  homogeneou*  equation* 

l  Ci  »  0 
j 


I  Cje  ■*  »  o 

j 

y  x2c.  -  r.  y  x.c.  =  o 
r  j  j  1  t  i  i 


f  (j  =  1,2,  3.  4) 


y  X j  Cj  e")  -  R2  y  XjCje^j  =  0 

j  j 


(2.3) 


Since  Xj  =  x,  Xz  =  - K,  X3  =  in,  X^  =  -ix  (eee  Appendix  I),  Equations  (2.3)  yield  the 
following  coefficient  determinant  which  must  vanish  in  order  for  these  equations  to  be  consistent. 


1 

e* 
x  -  R 


1 

(t  -  RJe" 


x  +  Rj 
(x  +  Rz)e-* 


.IX 


-(x  +  iRj) 
-(x  +  iRJe* 


-IX 


-<x  -  iRj) 

-(X  -  iR2)e-ilc 


=  0 


(2.4) 


Expanding  the  determinant,  we  have 

A  =  8i[2x2  sinh  x  sin  x  +  x(Rj  -  Rz)(sinh  x  cos  x  -  cosh  x  sin  x) 


+  RjRz(l  -  cosh  x  cos  x)  =  0 


(2.5) 


Given  the  appropriate  values  of  Rj  and  R2,  the  corresponding  values  of  x  may  be  determined, 
thereby  providing  the  necessary  values  for  =  Xj/L  =  x/L  for  the  frequency  equation 
[  Equation  (1. 1 3)] . 

The  related  displacement  equations  are 


K3KC  f  j/  xx\ 

Ujjjjjtx,  $,  t)  =  f  ^  cos  m$  sin  ^  2x  ysin  x  cosh  +  sinh  x  cos  J 
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+  kRj  ^r2  ^cosh  -  cos  +  I4  ^sin-y^  -  sinh  y^ 

+  2kR2  ^cos  k  sinh  y^  -  cosh  k  sin  yj  j  -  RjR2  £r3  ^o*1* 
-  cos  ^)  +  r2  (sin  £  -  sinh  yf)]J- 
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where 


2  . 


Tj  =  2k  sin  k  +  KRj(sinh  K  -  cosh  K  +  sin  K  -  cos  k)  +  2kR;>  cos  k 
+  Rj  Rgfsinh  k  -  cosh  k  +  sin  k  +  cos  k) 
r7  -  sin  k  +  sinh  k 


r3  =  cosh  k  -  cos  K 
r4  =  cosh  K  +  COS  K 


Having  defined  the  relation  for  the  deflections,  it  then  becomes  a  simple  matter  to  obtain  the 
accelerations  as  functions  of  time. 

For  the  case  where  the  elastic  restraints  at  either  end  of  the  cylinder  are  equal 
(R,  =  R2  =  R),  Equation  (2.  5)  reduces  to 


A  =  8i[x4  sinh  k  sin  k  +  R  (1  -  cosh  k  cos  k)]  =  0 


and  Equation  2.8  becomes: 


(2.7) 


(2.8) 


(2.9) 
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C  f  2  /  J(X  KX  \ 

wmn(x.+,t)  =  p-  cos  m$  sinc^nnt^K4  ^sin  it  sinh-£-  +  *inh  k  ein-£-J 


-  kR 

I"2 

/  KX 

\C04T  ‘ 

cosh^jj) 

+  r3 

(»inhT+  ^"t)] 

-  R2 

[r3! 

(cosh^ 

KX  \ 
-  co.rj 

+  r2  1 

(•i“T-  sinh^)]} 

where 


r'j  =  2k  sin  K  +  R(k  +  R)[sinh  K  -  cosh  K  +  sin  K  +  cos  k] 

Similar  expressions  for  umn(x,  t)  end  vmn(x>  t)  are  found  by  using  the  relations 


..  ..  K3  8wmn(x'<M> 

umn(X|**t)  =  Kj  5x 


,  ,  K2  8wmn<x-<M> 

vmn(x.4..t)  *--£■ 


T$~ 


(2.  10) 


(2.  11) 


If  the  cylinder's  end-closures  are  extremely  rigid  such  that  the  boundary  conditions 
may  be  assumed  fixed,  then  R  — oo  and  Equations  (2.  5)  and  (2.8)  reduce  to 


A  =  1  -  cosh  K  cos  K  =  0 


(212) 


wmn<x- 


,  .  C  f_  /  KX  .  Kx\ 

<M)  =  jmT  COS  m$  smu^t  j^r3  ^cob  -jj  -  cosh  -£•  J 

+  r2(sinh^  -  sin^)] 


(2.  13) 


where 


r"  =  sinh  K  -  cosh  K  +  sin  K  +  cos  K 

For  the  condition  where  the  cylinder's  end-closures  are  such  that  the  boundaries  are 
simply  supported,  then  R  =  0  and  Equations  (2.  5)  and  (2.8)  become 


A  -  sinh  K  sin  k  =  0  (2. 14) 

wmn  (x,  4,  t)  =  C4  cos  m$  sin  u^t  sin  ^  {Ref.  Eq.  (2.  3)J  (2. 15) 
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Equation*  (2. 12)  and  (2. 14)  provide  the  extreme  value*  of  k.  That  i*,  for  the  cylindri¬ 
cal  shell  with  both  edge*  fixed,  Equation  (2. 12)  give* 

K  =  j  [2n  +  1]  w  =  n'w  (n  :  1,2,3,,,.)  (2. 16) 


whereas  for  the  simply  supported  boundaries,  Equation  (2. 14)  gives 

k  =  nir  (n  =  1,2,3,...)  (2.17) 

Similarly,  from  Equation  (2.9),  we  obtain  the  variation  of  k  with  the  term  representing  the  end- 
closures'  stiffnesses.  This  is  shown  in  Figure  2. 1. 

2.1.2  Linear-Elastic,  Large  Deflection  Theory* 

In  order  to  avoid  the  prohibitive  difficulties  associated  with  a  detailed  analysis  and 
thereby  obtain  a  solution  that  will  lead  to  useful  information,  we  assume  that: 

(a)  The  elements  normal  to  the  middle  surface  remain  normal  and  unstretched 
(extensional  or  membrane  vibrations) 

(b)  The  cylinder  wall  is  thin  ^say  ^ 

(c)  The  longitudinal  (u)  and  tangential  (v)  displacement  components  are  small 
compared  to  the  radial  displacement  (w). 


For  the  equations  of  motion,  we  have 


8Nxy 
+  8y 


8N 

8y 


phU 


phv 


(2.18) 


Following  the  format  of  Fung  and  Sechler^'  we  define  the  strain  components  at  the  shell's 
middle  surface  as 


8x  + 

2  \  8x  / 

®v 

1  (izf 

w 

8y 

2  \  8y  /  ' 

a 

8u 

.  8v  .  8w 

8w 

8y 

+  8x  +  8x 

ey 

By  neglecting  the  inertia  terms  in  the  x  and  y  directions  (this  corresponds  to  a  first  order  approxi¬ 
mation  in  the  perturbation  procedure)  and  introducing  the  stress  function  F(x,  y)  where  Equa¬ 
tions  (2. 17)  are  satisfied  by 


8ZF 

8y2 


N 


xx’ 


N. 


yy* 


82f 

8x8y 


=  N. 


xy 


♦Notation  in  this  section  same  as  that  given  in  Appendix  I  except  where  noted. 
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FIGURE  2.  1.  VARIATION  OF  FREQUENCY  PARAMETER 
WITH  STIFFNESS  PARAMETER 
(See  Eq.  2. 9) 
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1 

we  obtain  the  compatibility  equation  1 


_4_ 

V  F 


1^ 

a 


and  the  equation  of  motion  for  the  radial  displacement 


4 

V  w 


■*£ 


82w 

7? 


a2„  .  2  2  .2 

8  F  8  w  8  F  8  w 

fcxBy  6x8y  +  0y2 


82F 

8x2 


1  8ZF 

+  I  77 


+  P<x.<fr. 


t)  -  pvfr  j 


(2.  20) 


(2.21) 


The  problem  now  becomes  one  of  assuming  a  form  for  the  radial  displacement,  using 
it  to  determine  the  stress  function  from  Equation  (2.20)  which,  in  turn,  is  applied  to  Equa¬ 
tion  (2.  21).  The  latter  equation  suggests  that  the  radial  displacement  may  be  of  the  form 

w(x.y,t)  =  f(t)  •  C(x,  y) 


where  G(x, y)  denotes  the  mode  shape. 


The  moat  appropriate  expression  is  one  which  is  sufficiently  complete  so  as  to  embrace 
all  the  significant  parameters  but  not  unduly  complex  so  as  to  complicate  the  analysis  and  intro¬ 
duce  refinements  that  are  not  commensurate  with  the  required  accuracy  of  the  final  numerical 
answers. 

As  an  example  of  a  simple  mode  shape,  we  take 

w(x,  y,  t)  =  f(t)[  sin  ax  cos  py  +  «l»  sin  ax  cos  yy]  (2.  22) 

where  a  *  ir/L,  (3  *  1  /a,  y  =  m/a,  and  m  denotes  the  number  of  circumferential  half  waves. 

Using  Equation  (2.  22)  in  Equation  (2.  20)  and  after  extensive  algebraic  manipulation,  we  arrive  at 
the  following  expression  for  the  stress  function 


F 


Eh  l  f(tfa2 
•  4 


[Kj  cos  2(Jy  +  K2  cos  2yy  +  K3  cos  (P  -  y)y  +  K4  cos  (P  +  y)y 


-I-  Kg  cos  2ax  4  Kg  cos  (p  -  y)f  cos  2ax  +  Kj  cos  (p  +  y)y  cos  2ax 

+  Kg  sin  ax  cos  Py  +  Kg  sin  ax  cos  yy)  (2,  23) 
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For  •  radial  displacement  equation  o f  the  form 

w(x,  y,  t)  *  f(t)[sin  ox  cos  Py  +  sin  ox  cos  yy  +  *  sin  ox] 
the  stress  function  is  as  follows 


F 


Ehffft)]  V 
4 


[Kj  cos  2py  +  K2  cos  2yy  +  Kj  cos  (P  -  y)y  +  cos  (p  +  y)y 


+  Kg  cos  2yx  +  Kg  cos  (P  -  y)y  cos  2ox  +  Ky  cos  (p  +  y)y  cos  2ax 
+  Kg  cos  2ax  cos  Py  +  Kg  cos  2ax  cos  yy  +  Kjq  cos  Py  +  Kjj  cos  yy 
+  Kj2  sin  ax  cos  py  +  Kjj  sin  ax  cos  yy  +  Kj^  sin  ax] 


where  Kj  through  K y,  Kj2  aod  Kjg  are  the  same  as  those  given  for  Equation  (2.  23)  (except 
that  4*  =  1 )  and 


Kg  = 


2*p 


..  2  .2.2 

(4a  +  p  ) 


K 


2* 

io=  pr 


Kg  = 


2*7 


(4a2  +  P2)2 


2* 

Kn  =  7* 


As  a  final  example,  for 

w(x,  $,  t)  »  f(t)(A  sin  ax  cos  Py  +  sin  ax  cos  yy  +  x  sin2  ax] 
the  stress  function  becomes 


F  =  — —  [Kj  cos  2Py  +  K2  cos  2yy  +  Kj  cos  (P  -  y)y  +  K^  cos  (P  +  y)y 

4  Kg  cos  2 ax  4  Kg  cos  (p  -  y)y  cos  2 ax  4  Ky  cos  (P  4  y)y  cos  2 ax 
4  Kg  sin  ax  cos  Py  4  Kg  sin  ax  cos  yy  4  KjQ  sin  3ax  cos  Py 
4  Kjj  sin  3ax  cos  yyj 


where  Kj  through  K^,  Kg,  and  Ky  are  identical  to  those  for  -Equation  (2.  23)  (for  >|»  =  1)  and 


k5  = 


T2  +  (PA)2-jfe 

8a4 


-[l-«y2Xf(t)I 
f(t)a(a2  4  y2) 


(224) 


(2.25) 


(2.  26) 


(2.27) 
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^  «t)a(az  +  pz)z 


jjsii 

(9a2  +  tT 


Kio  *  rr 


V”  =  (9a2  +  pV 


The  Gale r kin  procedure  may  now  be  used  to  find  the  differential  equation  for  the  time 
function  f(t),  thereby  providing  not  only  the  deflection,  but  also  the  acceleration  as  a  function  of 
time.  The  Galerkin  equation  will  have  the  form 


L  >»f  r 82w  82f  82w  82f 

/  /  w '  L  W  5^57 

+  p(x,  9,  t)  -  p#  J  w(x,  «f».  t)^  dxdy  =  0 


a2w  8ZF  2.  8Zf 

"  ~  8?"  a  8  x2 


8y 


(2.  28) 


where,  for  example,  the  stresa  functions  (F)  derived  from  the  radial  displacement  expression 
w(x,  t)  would  be  those  from  Equations  (2.  22)  and  (2.  23),  (2.  24)  and  (2.  25),  or  (2.  26)  and  (2.  27). 

In  order  to  obtain  a  solution  for  Equation  (2.  28),  it  is  first  necessary  to  specify  the 
external  loading  conditions,  p(x,  >)>,  t).  One  satisfactory  relation  is  the  following 


p(x,  «>,  t)  =  (p„  +  pj  cos  £  )  T(t) 


(2.29) 


where,  for  example, 
r(t)  =  e-6t 


or 


T(t)  =  ct 

The  differential  equation  for  the  time  function  is  of  the  form 

3  +  Aj  J3  +  A2J2  +  A3J  =  A4T(t)  (2.  30) 

wher£  J  ■  f(t)/a  and  A^, . .  .  are  coefficients  reflecting  the  particular  characteristics  of  the 
mode  shape  and  load  function.  As  a  particular  example,  we  take  L  =  36  inches,  A  =  6  inches, 
h  *  .036  inch,  E  =  2.9  X  10^  psi,  v  =  0.  25  and  p  =  24.7  X  10*^  lb-*ecZ/in^  as  the  properties  of 
the  cylindrical  shell.  For  the  mode  shape  defined  in  Equation  (2.26)  (with  x  =  0,  A=  1.0  or  0, 
and  n  *  2)  and  the  loading  defined  as 

p(x,  9,  t)  =  p0(l  +  B  cos  Py)e*6t 


the  numerical  vsdues  of  the  coefficients  are  as  follows: 


10 


For  A=  1:  p0  =  70  pel,  B  =  4/7  and  6  =  2,77  X  10*2 

Aj  =  23.4,  A2  *  -14.3,  Aj  =  5.8,  A4  =  3.72  X10’3 
For  A=  0:  p0  =  70  pai,  B  =  4/7,  and  6  =  2.77  X  10*2 

Aj  =  10.  4,  A2  =  -1. 1,  A3  =  0.  3,  A4  =  1.  65  X  10'3 


The  variation  of  J  aa  a  function  of  T(t)  ia  given  in  Figure  2.  2. 
diaplacementa  for  A=  1  and  A  =  0  are  reapectively 


0.022" 

0. 11" 


Accordingly,  the  maximum  radial 
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2.  2  Dynamic  Response  of  Circular  Plates 
2.2.1  Large  Amplitude » 


LIST  OF  SYMBOLS 


a 

radius  of  circular  plate  (in. ) 

D 

S 

12(1  -  v2) 

E 

modulus  of  elasticity  of  plate  (lb/ in2) 

f.  F,  G 

functions  of  time 

h 

thickness  of  plate  (in. ) 

•*n»  *n«  Yn, 

*n 

Bessel  functions  of  order  n 

N 

= 

(Nr  +  Nt)/(1  +  v)  (lb  in'1) 

N* 

buckling  load  of  circular  plate  (lb  in*1) 

Nr.  Nt 

stress  resultants  in  the  radial  and  tangential  directions,  respectively  (lb  in'1) 

P 

natural  frequency  of  plate  (sec*1) 

r.  e 

cylindrical  coordinates 

~r 

S 

£ 

a 

t 

time  (sec) 

T 

= 

th  /  E 

*2  V  12(1  -  v*)p* 

V 

velocity  (in.  sec*1) 

V 

ac 

y»2  /12(1  -  vz)p* 

h2  V  E 

w 

lateral  deflection  of  plate  (in. ) 

w 

* 

w 

h 

p 

mass  per  unit  area  of  plate  (lb  sec^  in' 3) 

p* 

mass  density  of  plate  material  (lb  sec^  in**)  • 

V4 

* 

2  2  2 

(V2)2  «  biharmonic  operator  *  (  -2—  +  —  ♦  -i- -2—  ^ 

\8r2  rTr  r2  8g2/ 

It  la  wall  known  that  whan  tha  lateral  deflections  ara  large  the  behavior  of  plates  is  gov* 
erned  by  two  coupled  nonlinear  partial  differential  equations  usually  credited  to  von  Kdrmdn.  (*•  2) 
The  dynamic  analogue  of  these  equations  has  been  derived  by  G.  Hermann.  (2.  3)  Chu  and 
Hermann^,  d)  have  calculated  the  fundamental  frequencies  of  rectangular  plates  at  large  amplitudes 
by  an  energy  procedure. 
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Howe  vs  r,  a  theoretically  exact  solution  of  the  dynamical  equations  even  for  the  case  of 
free  vibrations  is  nonexistent.  The  case  of  forced  vibrations  had  apparently  not  been  studied  at  all. 

It  is  to  be  realised  that,  owing  to  the  nonlinearity  of  the  equations,  the  free  and  forced  vibration 
problems  cannot  be  solved  separately  and  superposed,  neither  can  advantage  be  taken  of  the  concept 
of  normal  modes  as  is  usual  in  the  linear  case. 

In  1955,  Bergeit2,  5)  proposed  a  set  of  equations  which  may  be  termed  a  simplified  version 
of  the  von  Kdrmdn  system.  The  merit  of  his  method  lies  in  decoupling  the  two  equations,  so  that  one 
of  them  assumes  a  quasi- linear  form  and  can  readily  be  integrated.  He  showed  further  that  his  equa¬ 
tions  yielded  solutions  remarkably  close  to  those  obtained  by  more  elaborate  procedures  using  von 
Kdrmdn's  equations,  at  least  as  regards  deflections.  Berger,  it  must  be  remarked,  confined  him¬ 
self  o  the  static  case. 

Nash  and  Modeert2*  extended  the  Berger  equations  to  the  dynamic  case  for  rectangular 
plates  and  showed  that  the  free  vibration  problem  can  be  treated  to  yield  results  which  are  reason¬ 
ably  close  to  the  results  obtained  by  Chu  and  Hermann.  (2*4) 

In  what  follows,  the  dynamical  equations  k  la  Berger  are  extended  to  circular  plates.  By 
using  an  approximate  step-by-step  procedure, it  is  shown  that  not  only  the  free  vibration  but  also  the 
forced  vibration  problem  may  be  solved  using  the  elementary  concept  of  normal  modes.  It  is  neces¬ 
sary  to  add  that  although  the  basic  Berger  formulation  has  a  quasi-linear  form  for  the  static  case,  it 
is  nonlinear  for  the  dynamic  case  and  cannot  be  solved  with  comparable  ease. 


2.2. 1.1  Equations  of  Motion 

For  the  dynamical  case,  the  Berger  equations  for  circular  plates  may  be  written  in  the 

form 


V*w  -  a-Vw  +£-  i-i,  iiltitii 

D  gt2  D 


a2h2  flu  ,  1  /  8w\ 2  u,  1  Bv,  1  /  ®w\ 2  . 

—  -K+1KT7J  +  i*tT0+2^2  \Te)  "  +  et 


(2.31) 


in  which  V*  =  (V2)2  is  the  biharmonic  operator,  p  is  the  mass  of  the  plate  per  unit  of  area,  w  is  the 
lateral  deflection,  u  and  v  the  radial  and  tangential  displacements,  D  % 


Eh 


,  is  the  flexural 

12(1  -  v2) 

rigidity  of  the  plate,  h  its  thickness,  E  is  the  elastic  modulus,  v  is  Poisson's  ratio  and  f  a  given 
forcing  function. 


The  radial  and  tangential  strains  er  and  eg  may  be  expressed  in  terms  of  the  stress  result¬ 
ants  Nr  and  N*  as  follows 


er-^(Nr-vNt) 


ct^  (Nt-  vNr) 


(2.32) 


It  follows  from  (2.  32)  and  the  second  of  (2.  31)  that 


q2hZ 

12 


1  -  v  N  .  Nh‘ 
~~hE~  (Nr  +  Nt>  *  JJ5 


(2  33) 


In  the  Berger  formulation,  a  is  constant  throughout  the  plate  and  thus  N,  which  has  the 
same  dimension  as  a  stress  resultant,  is  also  constant  throughout  the  plate.  (The  simplicity  of  the 
Berger  equations  apparently  derives  from  the  basic  assumption  that  (Nr  +  Nt)  is  constant.  ] 
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From  Equation!  (2.  33)  and  the  aecond  of  (2.  31)  one  finde 
Nh2  8u  1  /  8w\2  u  ,  1  8v  ,  1  / 8w\  ^ 

=  07  +  I  1*7/  +  7  +  7  *e  +  -£i  [FeJ 


12D 


In  view  of  the  assumed  constancy  of  N  with  reepect  to  r  and  6,  one  may  multiply  by  r  dr  d0  and  inte¬ 
grate  over  the  plate  to  find 


is rdr de  =  //(It  *  t)  *  ffi  si  r,1"‘9 

4 XT  [(if)2  '-H  (li)2] 


If  the  boundary  conditions  are  such  that  u  and  v  vanish  at  the  boundaries,  the  first  two  integrals  on  the 
right-hand  side  are  evanescent  and  there  must  follow 


2ir  a 


Nh2ira2  1  f  C  /8w\2  ,  1  /  8w\ ^ 

-1Id-=  2  J  J  Is?)  +7z  \Se)  rdrde 


0  0 


using  (2.33)  the  first  of  (2.31)  may  be  written 

v4w-£v2w+|-  iiliAil 

D  D  8t2  D 

The  problem  now  reduces  to  the  integration  of  the  system  (2.  34). 


(2. 34) 


2.  2. 1.2  Integration  Procedure 

It  is  seen  from  the  first  of  (2.  34)  that  N  is  a  function  of  time  alone  and  that  the  second  of 
(2.34)  is  effectively  nonlinear.  The  problem  of  finding  the  dynamic  response  of  such  a  system  cannot 
be  solved  by  superposition  of  the  separate  solutions  of  free  vibrations  and  of  the  steady  state  prob¬ 
lem. 


A  step-by-step  integration  procedure  is  here  proposed  for  the  integration  of  the  system 
(2.  34).  The  assumption  is  made  that  N,  which  is  proportional  to  the  sum  of  the  membrane  tensions 
(Nr  +  Nt),  remains  constant  at  N  for  a  sufficiently  short  interval  of  time  A  .  It  is  evident,  that  the 
second  of  (2.  34)  assumes  a  linear  form  and  may  be  readily  integrated.  N  is  evaluated  from  the  first 
of  (2.  34)  at  the  beginning  of  each  time  interval.  The  nonlinear  problem  is  thus  reduced  to  a  series  of 
linear  problems  in  time. 

Even  this  simplification,  it  may  be  noted,  still  involves  a  prodigious  amount  of  labor.  For, 
since  the  normal  modes  of  vibration  are  a  function  N,  they  are  a  function  of  time  also.  Furthermore, 
the  forcing  function  f(r,  6,  t)  has  to  be  expanded  in  a  series  of  the  normal  modes  appropriate  to  that 
time  interval. 

Thus,  the  general  procedure  is  as  follows:  Starting  with  the  given  initial  conditions,  eval¬ 
uate  71  from  the  first  of  (2.  34).  Use  this  value  of  N  to  solve  the  homogeneous  part  of  the  second  of 
(2.  34)  and  determine  the  natural  frequencies.  Expand  the  forcing  function  in  terms  of  the  normal 
modes  and  determine  a  particular  solution  of  the  second  of  (2.34).  Use  the  initial  conditions  at  time 
zero  to  determine  the  deflection  and  velocity  at  the  end  of  the  time  interval  A.  Use  the  latter  as  the 
initial  condition  for  the  second  time  interval  A  Repeat  this  procedure  up  to  any  desired  time  t  =  sA. 

(Some  of  the  tediousness  of  this  procedure  may  be  reduced  if  one  ignores  the  change  in  the 
normal  functions  at  each  time  interval,  using  throughout  instead  the  functions  appropriate  to  the  initial 
interval,  but  faking  account  of  the  change  in  the  natural  frequencies  at  each  interval. ) 


15 


Returning  now  to  the  second  of  (2.  34),  let 
w  «  R(r)0(0)  tin  pt 


(2.35) 


where  R  is  e  function  of  r  alone  and  9  ia  a  function  of  6  alone.  On  substitution  into  the  homogeneous 
part  of  the  second  of  (2.34)  and  dividing  out  the  common  factor  ainpt  there  follows 


V4(R0)-  ^  V2(R0)  -  sh.  R0  *  0 
which  may  be  written  in  the  form 


where 


(R0) = 0 

(2.36) 

•2-¥  [(>^) 

1/2  "I 

-J 

► 

(2.37) 

1/2  -l 

“J  J 

»2  .2  ,  Na2 

P  -a 

(2.  38) 

(2.39) 

Equation  (2.36)  shows  that  the  complete  solution  may  be  obtained  by  adding  together,  with 
appropriate  arbitrary  constants,  the  solutions  of  the  two  equations 


and 


M) 


,2. 


R0  =  0 


RO  -  0 


(2.40) 


a  2/ 

On  letting 
*  0  =  coa(n6  -  7n) 

Equations  (2.40)  yield  two  ordinary  differential  equations  for  the  determination  of  R 

i!5+IdR+  (zl.  Si)  R  =  0 

dr2  r  dr  \  a2  r2/ 

+  ii*.  f£+  si)  R=0 

dr2  r  dr  \a2  r2  / 

The  solutions  of  (2.42)  are,  respectively 

R ■"*(¥)  *cr"(¥) 


(2.41) 


(2.42) 


and 
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RaBIn(¥)  +  DK»(¥) 


in  which  Jn  and  Yn  denote  Beeeel  function*  of  the  order  n  of  the  firet  and  eecond  kinds,  respectively. 
Iq  and  Kn  are  the  modified  Bessel  functions  of  the  first  and  second  kinds,  and  A,  B,  C,  and  0  are 
arbitrary  constants.  (The  notation  for  the  Bessel  functions  is  that  used  in  Reference  2.  7. ) 

In  the  case  of  a  solid  plate,  Yn  and  Kn  are  inadmissible  because  of  the  singularity  at  r  *  0, 
and  the  general  solution  of  the  system  (2.42)  is  thus 


R  =  A„Jn  1  BnIn  (2.43) 

The  general  solution  of  the  homogeneous  part  of  the  second  of  (2.34)  may  now  be  written 
"{-(?)*-■(?)}  {  cosn0  +  Cnsinnd|.  ^sinpnt  +Dn  cos  Pnt^  (2.44) 

The  frequency  equation  and  the  ratio  Bn/An  may  be  obtained  from  the  boundary  conditions. 


In  the  case  of  a  simply  supported  plate,  one  has 


w  =  0 

8^w  /  l_  8w  _1_  8  Zw\ 
8r2  V*  8r  r2  8  02  / 

and  in  the  case  of  a  clamped  plate 


=  0 


at  r  =  a 


w  =  0 
8w 


57 


=  o 


at  r  =  a 


(2.45) 


(2.46) 


Substitution  of  (2.44)  into  (2.45)  yields 
Bn  Jn(tt) 

An  =  '  WP) 

Jn  +  l(°)  In  +  l(P)  P2  +  u2 
°  Jn(<0  +  P  In(P)  1-v  . 


(2. 47)  ' 


and  the  substitution  of  (2. 


An  ‘  In(P> 


Jn  +  1*0> 
a  Jn^) 


+  P 


44)  into  (2.  46)  leads  to 


In  +  l<P> 
WP) 


(2.48) 


The  second  of  (2.47)  is  the  frequency  equation  for  the  simply  supported  plate  and  the  secofld  of  (2. 48) 
for  the  clamped  plate.  These  equations,  together  with  (2.3*8),  are  sufficient  for  the  determination  of 
the  natural  frequencies  Pn  for  given  values  of  TT. 
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Defining 


2L»  4 

N*  9 


N*  =  with  v  *  0.  3,  for  a  limply  supported  plate 


14.68D 


for  a  clamped  plate 


(2.49) 


the  frequency  equation  for  the  simply  supported  plate  may  be  written 
Jn  +  l<a>  ,  .  ^  +  1<M  P2  +  a2 

“  Jn(“)  1  *  v 

p2  -  o2  =  4.  24 
w  =  0.3 


(2.  50) 


and  that  for  the  clamped  plate 

Jn  +  1<°>  .  .  ^  4  1<P)  _ 

Jn(“)  P  In(P)  *  ° 

p2  -  o2  =  14.684 


(2.51) 


Equations  (2.  50)  and  (2.  51)  may  be  solved  for  various  values  of  the  parameter  4  and  the 
values  of  a  and  p  so  derived  may  be  tabulated  or  graphed  as  functions  of  4. 


B, 


hi  w 


It  will  be  noted  from  (2. 45)  and  (2. 46)  that  the  ratio  =  •“•••  J.  whether  the  plate  is 

Anj  In(Pj) 

simply  supported  or  clamped. 


Letting 


(2.52) 


one  may  write  the  deflection  function  (2.44),  corresponding  to  the  (nj )th  frequency  as  follows 


'*hj  =  Rnj  {AnjC°»nS  +  Bnjsinne}  {»inpnjt  +  Dnjcospnjt]- 


(2.53) 


in  which  Anj,  Bnj  and  Dnj  are  new  arbitrary  constants,  to  be  determined  from  the  initial  conditions. 

(2-54) 


mu  “nj  • 

It  may  be  seen  that  (2.  53)  may  be  written  in  the  form 
wnj  =  Rnj  { A„j(t)  cos  nfl  +  Bnj(t)  sinnS } 


where 


Anj(t)  ■  Anj  {sinpnjt  +  D„j  cos  Pnjt} 
Bnj(t)  ■  B„j  ^sinpnjt  +  Dnj  cos  Pnjtj 


(2.55) 


Any  arbitrary  deflection  satisfying  the  boundary  conditions,  and  therefore  the  solution  of  the  homo¬ 
geneous  part  of  the  first  of  (2.34)  appears  finally  in  the  form 

wc  =  £  £  RnjAni(t)cosn8  +  £  £  RnjBnj(t)  sinn  (2.56) 

n=  0  1  n= 1  j=  1 

where  the  superscript  "c"  on  the  left  denotes  "complementary  function.  " 
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2.  2.  1 . 3  Forced  Vibration 


In  »eeking  a  particular  solution  of  the  first  of  (2.  34),  it  will  be  supposed  that  the  forcing 
function  f(r,  0,  t)  can  be  expanded  in  terms  of  the  principal  modes  in  the  form 


f(r,  0,  t)  =  F(t)Z£Rnj(rnj  cos  n0  +  snjsinn0) 
A  particular  solution  of  (2.  34)  is  sought  in  the  form 
wP  s  SSRnj  ^anj(t)  cos  nfl  +  bnj(t)  sinnsj 
in  which  anj  and  bjjj  are  undetermined  time  functions. 


(2.57) 


(2.58) 


Inasmuch  as  Rqj  cosnS  and  R^j  sin  n0  satisfy  the  homogeneous  part  of  the  first  of  (2.  34) 
there  must  follow 


/  cosn0\  /  cos  n8  \  2  /  cosn0\ 

V4  (Rnj  -  ?V2  Kj  )=^  Kj  ) 

\  sinnd  /  D  \  sinnfl  /  ®  \  sinnO  / 

Substituting  (2.  58)  into  the  first  of  (2.  34)  and  making  use  of  (2.  59)  one  finds 


(2.59) 


..  ,  2  FM*nj 

anj  +  Pnj'nj  " 


P 

F(t)i 
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(2.60) 


bnj  +  Pnjbnj  =  p 

Let  the  particular  solutions  of  these  equations  be 

anj  =  Unj0*4* 
bnj  =  bnjG(t) 


} 


(2.61) 


From  (2.  58)  and  (2.61),  one  gets 

wP  =  (gn:  cos  n0  +  hnj  sinn0)G(t)  (2. 62) 

n  j  J  J  J 

The  complete  solution  of  the  first  of  (2.  34)  may  now  be  obtained  by  adding  (2.  56)  and  (2. 62) 

co  co  CD  oo 

w  =  wc  +  wP  =  £ 


n^Q  R^jlAnjW  +  gnjG(t)]  cos n0  +  ^  +hnjG(t)l  sinn0 

(2.63) 


Let  the  initial  conditions  be 


w  =  4><r)+<0)  "I 

*  =  x(rMe)  J 


at  t  =  0 


These  can  always  be  expanded  in  Fourier-Bessel  series  of  the  form 

GO  GO  GO  OO 

w=  £  £  RojHnjcoen0+  £^  £^  RnjLnjsinn0 


oo  co 
w  =  £ 


1  n  =  0  Rn'iMn'i  C°*  +  j  ?  1  n  =  0  Rn-’N"'i 


sinn0 


(2.64) 
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where 


$(r  H^flJRoj  coa  nfl  rdrdfl 


^(r)4i(0)Rnj  sin  nfl  rdrdfl 


x(r)p(0)Rnj  cos  nfl  rdrdfl 


n  i  0 


*Nnj  = 


rr 


X<r)fi(9)Rnj  sin  nfl  rdrdfl  n^O 


(2.65) 


in  which 


X  = 


R„jr  dr 


n 


=  0 


=  w/  R 

0 


nj 


r  dr 


n  i  0 


From  (2.  63)  and  (2.  64),  using  (2.  65),  one  finds 

QO  CO  CO  CD 

w=  2  2  RjjjU,,;  cos  nfl  +  2  2  Rn:Wnjsinn6 

n  =  0  j  =  1  n  =  1  j  =  1 

oo  oo  oo  oo 

w=  2  2  RnjVnj  cos  nfl  +  2  2  RnjZnjsinnfl 

n=0  j= 1  J  n  =  1 j  =  1 


(2.66) 


in  which 


Unj  =~  {Mnj  ’  «nj6<°)}  •‘"Pnj*  +  { «n j  *  *njG<0>}  co»  Pnj‘  +  8njG<‘> 

Vnj  =  {Mnj  -  8nj6<°>}  co,Pnj‘  *  Pnj  {Hnj  *  8njG(0)}  sinpnjt  +  gnj6(t) 

Wnj  “  {Nnj  *  hnjG<°>}  "‘"Pnj*  +  {Lnj  ‘  hnjG<°>  }  cos  p„.t  +  hnjG(t) 
nj 

Znj  =  {Nnj  *  hnjG<°>}  coi  Pnj‘  '  Pnj  {Lnj  "  hnjG(0>  }  *in  Pnj*  +  hnjG(*> 

The  deflection  is  now  completely  known  and  the  acceleration  may  also  be  determined.  It 
is  to  be  remembered,  however,  that  this  solution  is  valid  only  for  the  time  integral  (t  -  r)  during 
which  the  parameter  R  is  supposed  to  remain  constant.  In  Equation  (2.66),  therefore,  the  time  vari¬ 
able  must  be  (t  -  t)  instead  of  t.  The  parameter  N  is  determined  by  the  first  of  Equations  (2.  34)  with 
t  taken  at  the  beginning  of  the  time  interval  (t  -  t).  Thus 


6D 


N  = 


wa2h2 


2ir  a 

/  /  [(S?V 


(If)2] 


rdrdfl 


t  =  r 


(2.67) 
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Using  (2.  63)  and  (2.  66)  and  (2.  50)  and  2.  51),  Equation  (2.  67)  may  be  written 


■rdr  + 


[2/  {jf,  J, 

L  0 

+  /  {  SEE  RijIMxktHnjHnk  +  LnjLnk))  rdr  + 
q  l  n= 1  j  = 1  k=  1  j 

+  f  {  n?1  kfj  n^RnjRnk(HnjHnic  +  Lnjl^k) }  ^ 

ft  J 


in  which  R^j 


(2.68) 


The  integrals  in  (2.68)  cannot,  unfortunately,  be  readily  obtained  in  closed  form,  and 
recourse  must  be  had  to  numerical  integration. 

2.2.  1.4  Summary  of  Integration  Procedure 

In  numerical  integration  it  is  advantageous  to  use  nondimensional  quantities.  These  will  be 
distinguished  by  placing  a  bar  over  the  corresponding  symbols.  Let 


"t  =  ±  /  n::z 

aZ  V  12(1  -  v2)p* 

where  p *  is  the  mass  density  of  material  of  the  plate 


£ 

a 


Let  s  denote  the  number  of  intervals  A  from  time  zero.  Then  the  deflections  and  velocities  at  the  end 
of  the  s™  interval  are  given  by 


*s 


CO  00  OO  CD 

2  2  R_;  Uni  cos  n©  +2  2  R»,j  Wn:  sinnO 

j=ln  =  0  J>  J'  j=ln=l  J»  •’» 


00  00 
^8=2  2 
j  =  1  n  =  0 


oo  oo 


>  s  >  i 


j 


(2.69) 


The  Rnj  are  defined  by  Equation  (2.  52)  and: 
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Unj.  *  p^-  {“*.  '  *nJ.6tA(‘  •  l)]}  ,in(pnjBA)  + 

+  {^j,  *  *  Ml}  co.^A)  +  g^OlAe] 

Vnj,  =  {^j.  *  *nj,6Ut*  *  Ml}  co.(pnj>A)  - 

'  P»j.  {  “aj,  ‘  *nj,G[AU  '  Ml}  •*a<Pnj,A>  +  8nj§Gt  A,1 

Wnj.  ■  Ki.  ‘  A(*  *  »«}  *>***'  + 

+  {^j,  *  NijgG^  -  l>l}  cos(pnj#A)  +  h^sG[A«] 

Znj„  =  {Nnj„  *  hnjgGlA<a  *  Ml}  co.  (pnj^A)  - 

•  P“j.  {Si.  -  hnjBGtA(*  *  1)]}  ,in(pnj.A)  + 

* 

and 

Hnj,  =  £  j  £  Unk>  x  Rnk^  _  j  Rnj>rdr  j £  R^rdr,  n  >  0.  s  >  i 

Mnj.  =  Jj  £  ^nkf  _  j  ®nk#  .  j  Rnj,7d?  jj  R^j8Tdr,  n  >  0.  .  >  i 

=  k*  j  /  Wnj,  .,*Bk|.1  *nj  .7dF//  R^.7dF>  n  >0>  *  >  : 

Nnj<  =  %x  £  Ynk,  _  j  R„k>  _  j  RnjJdT  R^Tdr,  n  >  0,  .  >  2 

and  Hnjj,  L,^,  Mnjj ,  and  Nnj  ^  are  defined  by  the  initial  condition*. 

The  natural  irequencie*  are  given  by 


a>l  (2.70) 


(2.71) 


pnj#  =  Pnj,  anj. 


(2.72) 
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p.2  -  • l  -  12  J  ,  R4j.  ^k.  Hoj.  HokJ  rdr  + 

+  6  /  {  Z  Z  a  Rfcj  Rnk.  (Hnj.Hnk,  +  Lnj.Lnk,)}  ?dT  + 
0  L  n  =  1  j  =  1  k=  1  J 

+  6  /l  {  ji,  j?,  k?!  n2Rnj.  *■*.  (H"is  «■*.  +  w}  T 


+  ,  s  >  1 

(2.73) 


°JS  d? 

Some  of  the  above  relations  may  be  simplified  by  noting  that  if  A  is  taken  sufficiently  small 
Rnjg  _  j  *  Rnj 

I  Rnkg  _  !  Rnjg  »«  *  °  k  ^  j 
Equations  (2.  71)  then  become 


^8  “  U-J.  -  1 

n  >  0 

M«VjB  "  Vnjg  .  ! 

n  >0 

^8  "  WnJ8  -  1 

n  >  0 

Nnjg  *  Y"j8  .  ! 

n  >  0 

Equation  (2.73)  simplifies  to 

P2  -  a2  =  12  SZ|ojkHoj9  Hok(j  +  6  SZZ|njk(Hnjs  1*^) 

+  6  ZZZri..  n2(H  .  H.  +  L  .  L,  )‘ 
nj  k  ''jk  nj8  J'jg  nk8 


(2.  74) 


(2.75) 


^ojk  =  /  R0jo  Rok0  rdr 
^njk  =  /  Rnj.  Rnk-  rdr 


Vik  =  /  Rnkn  r 


o  “•’o  o  r  J 

If  the  right-hand  side  of  (2.  75)  is  evaluated  for  any  value  of  s,  then  setting 
2  2 

PJ  -  a,  =  4.  2  <(>  for  simply  supported  plate 

=  14.  68  4  for  clamped  plate 


(2.76) 
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on*  may  determine  4  and  thereby  tha  valuaa  of  aBj  and  pnj  If  these  have  been  tabulated  for  varioua 
values  of  +•  (See  note  under  Equation  (2.  51 ).  ]  Equation  (2.  72)  may  then  be  uaed  to  determine  the 
nondimenslonal  frequencies. 

2.  2. 1.5  Numerical  Example 

The  response  of  a  circular  plate  to  a  loading  pul  so  will  be  investigated  by  the  method 
described  on  the  preceding  page. 

A  blast  of  short  duration  is  best  treated  as  an  impulse.  Assuming  that  the  blast  is  deliv¬ 
ered  uniformly  over  the  area  of  the  plate,  let 


/( 

f(t)  dt  =  I 

e  being  the  time  at  which  the  pulse  decays  to  zero. 

By  the  momentum  principle,  the  uniform  starting  velocity  is 


(2.77) 


where  M  is  the  total  mass  of  the  plate. 


The  plate  under  consideration  has  the  following  geometric  and  material  properties 
a  =  6",  E  =  30  X  106  psi,  p*  =  7.  34  X  10'4  lb  sec2  in'4 


Thus, 


- 

h2  V 


!l£i 


*  0. 45  sec  in' 


The  initial  velocity  v0,  calculated  from  the  characteristics  of  the  charge  (used  in  an  experiment,  the 
details  of  which  are  omitted  here)  and  relation  (2.  77)  above  gives  an  estimated  initial  velocity  of 
893  in.  sec'*.  Thus,  the  nondimenslonal  initial  velocity  is 

vD  =  893  X  0.45  M  400 


The  response  of  the  plate  will  be  analysed  on  the  assumption  that  the  plate  receives  a  uniform  velocity 
of  400  while  at  rest. 


The  parameter 


_h_ 

.2 

Thus 

t  =  61. 6  t 

This  gives  the  relation  between  nondimenslonal  time  T  and  true  time  t.  The  nondimenslonal  time 
interval  A  will  be  taken  as  0.0616.  This  corresponds  to  a  true  time  interval  of  one  millisecond. 

The  basic  equations  given  in  Section  2.  2.  1.4  can  be  simplified  considerably  because  of 
axial  symmetry.  In  what  follows,  the  bars  over  the  nondimenslonal  quantities  will  be  dropped.  No 
confusion  can  arise  as  all  qusmtities  henceforth  are  nondimensional. 


E 


12(1  -  v2)p* 


=  61.6 
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Equation*  (2. 69)  become 
**  =  j?i  *^1  U°j» 

=  j?!  ^Ji  Voj. 


>  1 


(2.  78) 


From  Equation*  (2.  70),  (2.  71)  and  (2.  74) 

Hoj8  "  Uoj8  .  !  Bin  (poj8  .  ,A>  +  Hojs  .  ,  co,<Poj8  .  A> 

Ja  -  1 

Mojs  =  Voj8  .  ,  =  Moj„  .  ,  “« <PojB  .  !  A>  -  Poj8  .  !  Hoj8  .  !  sin<Poj8  .  J 
The  initial  conditions  give 
Hojj  =0 


A) 


*  >  2 


(2.79) 


M°j1  =  v° 


/  Roj  rdr  /f  r|j 

0  /  0  1 


rdr 


where 


with 


Thu# 


Jo(°j  ) 

Roj  =  Jo(aj,  r)  -  T_/„.  i  *o(aj,  r> 


ajj  =  2. 222  for  j  =  1 


=  5. 452  for  j  =  2 


P  =  «?, 


%=i~8in(P°jA> 


r°Jl 


Voji  =  MOjj  co»  <PojA> 


From  Equation  (2.  72) 

p  .  =  a  .  B  . 
°J8  °J8  °J8 


s  >  1 


(2.80) 


(2.81) 


Using  j  =  1,  2  only,  Equation  (2.  73)  become* 

U  fe) 

(tt? 

Assuming  that  the  plate  is  simply  supported 
p2-a2=4.2*8 


1  dR0i  dR 


rdr  +  24 


r  h  h 

H  H  J  ——1  -  ■ --* 

°1  °2,  »  dr  dr 


rdr  + 


rdr  s  >  1 


(2.82) 


Thus,  68  known,  and  a(  and  P8  may  be  read  from  chart*  already  prepared. 
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2.2.2  Plastic  Collapse 


LIST  OF  SYMBOLS 


v 

U 

R 

• 

R 

a 

h 

P 

t 

r 


'o 


H 

ffr 

<rfl 


er*  '6 
Z 

v 

E 

P 

P 

n 

c2 

6 

K 


initial  velocity  given  to  membrane,  a  constant 

radial  velocity  of  flat  portion  of  membrane,  positive  outward 

radius  of  flit  portion  of  membrane 

radial  velocity  of  bending  wave 

initial  radius  of  membrane 

initial  thickness  of  membrane 

mass  density  of  membrane  material 

time 

radius  of  position  of  any  particle  in  the  flat  portion  of  membrane  at  any  time 

initial  radius  corresponding  to  r 

thickness  of  membrane  at  any  time 

radial  stress 

circumferential  stress 

yield  stress  in  teneion 

natural  (logarithmic)  strains  in  the  radial  and  circumferential  directions 
distance  of  central  flat  portion  from  initial  plane  of  the  diaphragm 
Poisson's  ratio 
modulus  of  elasticity 
=  R/a 
=  U/v 
=  H/b 
=  <r/p 

=  vt/a  =  Z/a 
=  4<r/pv2 

=  [(4<r/pv2Ml  -  v)/E]  */2 


In  investigating  this  problem  it  was  discovered  that  there  are  certain  existing  solutions  to 
corresponding  problems  in  underwater  shock  which  have  been  satisfactorily  verified  experimentally. 
While  the  problems  of  blast  effects  in  air  or  of  soft  radiation  are  by  no  means  the  same  as  the  effects 
of  blast  loads  underwater,  they  are  sufficiently  similar  to  warrant  the  supposition  that,  as  a  first 
approximation,  the  two  problems  may  be  studied  by  identical  procedures. 
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The  eolution  to  the  pleetic  collepee  of  a  diaphragm  proposed  by  Hudson(2-  8),  and  a  similar 
solution  by  Frederick^'  9),  refer  specifically  to  underwater  shock.  It  is  proposed  now  to  use  the 
results  of  these  analyses,  particularly  that  by  Hudson,  to  predict  the  collapse  of  plates  due  to  blast 
effects  in  air.  These  predictions  will  then  be  subjected  to  experimental  checks  under  the  test  pro¬ 
gram  at  present  being  formulated.  If  it  is  found  that  the  discrepancies  between  theoretical  and 
experimental  results  are  sufficiently  serious,  a  fresh  theoretical  analysis  will  then  be  undertaken.  It 
is  believed  that  this  approach  will  lead  to  the  greatest  economy  of  effort  in  the  overall  program. 

2.  2.  2.1  Summary  of  the  Hudson  Approach 

The  problem  attacked  by  Hudson  is  the  damage  done  by  underwater  explosions  to  thin  metal 
circular  diaphragms,  air  backed,  and  held  rigidly  at  the  peripheries.  The  material  of  the  diaphragm 
is  supposed  to  be  rigid-perfectly  plastic.  Although  Hudson  has  also  given  an  approximate  solution  for 
a  work-hardening  material,  the  primary  difference  noted  is  that  in  the  work-hardening  solution  the 
apex  of  the  cone  {the  final  shape  assumed  by  the  diaphragm)  is  rounded  off,  whereas  in  the  perfectly 
plastic  material  the  apex  is  a  point.  For  the  present  purposes,  it  appears  sufficient  to  fix  attention 
on  the  solution  for  the  rigid-perfectly  plastic  material. 

The  shape  assumed  by  the  diaphragm  at  any  time  t  <  t8,  where  t8  is  the  total  time  for 
deflection  or  "swing  time,  "  is  that  of  a  truncated  cone  as  shown  in  Figure  2.3.  The  flat  central 
region  travels  with  a  constant  speed  v  normal  to  its  plane,  decreasing  in  radius  with  its  motion,  the 
diaphragm  finally  assuming  the  shape  of  a  cone  with  a  zero  thickness  at  the  apex. 

While  an  explicit  form  of  the  solution  of  the  equations  of  motion  is  not  given  even  for  the 
highly  idealized  model  proposed  by  Hudson,  certain  special  cases  can  be  solved  and  information 
gleaned  regarding  the  following: 

(1)  The  radius  R  of  the  bending  wave  (see  Fig.  2.  3)  as  a  function  of  time 

(2)  The  diaphragm  profile  at  each  instant 

(3)  The  thickness  distribution 

(4)  Displacement-time  curves  of  particles  in  diaphragm 

(5)  Stress  and  strain  distributions 

(6)  The  center  deflection  as  a  function  of  v 

(7)  The  total  time  for  the  deformation  to  take  place 
The  solutions  given  by  Hudson  are  summarized  below. 

2.  2.  2.  2  An  Elementary  Approximation  to  the  Solution 

When  the  acceleration  of  the  material  in  the  flat  central  portion  is  negligibly  small  through¬ 
out  the  motion,  i.e. ,  U  is  a  constant,  the  solution  is  as  follows 


(2.83) 


U  =  v2/2c 

(a) 

R  =  a  -  ct 

(b) 

/  r  \  (v/c)2 

.  (c) 

Z  =  -  (a  -  R) 

(d) 

The  thickness  distribution  in  the  deformed  diaphragm  given  by  (2.83  c)  shows  a  dimpling 
tendency  at  the  center.  In  fact,  at  the  last  moment,  the  thickness  becomes  zero  at  the  center.  How¬ 
ever,  as  the  variation  of  the  thickness  with  R  is  so  rapid  near  R  =  0,  the  tiny  pinhole  may  not  be  very 
apparent  experimentally. 
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FIGURE  2.  3.  A  SKETCH  OF  THE  DIAPHRAGM  CONFIGURATIONS 
AT  TWO  SUCCESSIVE  INSTANTS,  t  AND  t+dt.  THE  VARIOUS 
QUANTITIES  WHICH  ENTER  INTO  THE  SEVERAL 
GEOMETRIC  AND  KINEMATIC  RELATIONS 
ARE  PORTRAYED  HERE 
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Equation  (2.  83  d)  shows  that  the  diaphragm  auumei  a  conical  shape  whose  center  deflec¬ 
tion  ie  proportional  to  the  initial  velocity  v. 

The  total  time  of  deflection  t(  is  given  by 

t.  =  a/c  (2.84) 

2.  2.  2. 3  The  Exact  Solution  for  a  Rigid- Perfectly  Plastic  Material 


where  xj. 


i(xi  -  l)/(xi- x2)  -|(x2- l)/(x2- xj) 


R 

/  H2  -  x 1 

L  9  ^ 

\  ^  /  H-2  -  X1  \ 

(a) 

7  ■  PM  = 

(  X2  -  Xj 

)  (x2-x2) 

-2(x1)/(x1  -  x2)  ,  -2{x2)/(x2  -  xj  ) 

r 1  *■*> =  1 

(  ^2-xl 

\  X2  -  Xj 

)  *(£*)  ' 

(b) 

f  =  4(p)  = 

4 

/  (n2- 

-bL. — •  P(p)dp 

X1 ) (p2  -  *2) 

(c) 

X 

x2  are  the  smaller  and  larger  (positive)  roots  respectively  of  the  quadratic  equation 
x2  -  (26  -  3)  x  +  2  =  0  (2.86) 

It  is  to  be  noted  that  the  initial  conditions  are 


p=t|=l,  p  =  x  when  4=0 

Equation  (2.  85  c)  is  not  integrable  in  terms  of  elementary  functions. 


3  3  /■“ 

occur  in  pairs  of  positive  values  for  —  +  V2  <6  <. oo  and  as  pairs  of  negative  values  forO  <^6  <—  -  V2 


For  positive  values  of  6,  which  are  the  only  ones  of  physical  interest,  the  roots  of  (2.  86) 

3 
•2 

(a  very  small  range).  For  intermediate  values  of  6,  the  roots  are  complex.  Hudson  confines  himself 
to  a  consideration  of  values  of  6  2-^  +  \  Z  .  At  the  lower  limit  of  this  range,  corresponding  to  very 

large  v,  or  small  tr,  values,  we  have  xj  =  x2  =  ‘Jz  .  As  6  becomes  very  large,  'll  >  xj-~  0  and 
«/2  <  x2  -»  ao. 

A  case  of  considerable  interest  arises  when  k2  coincides  with  one  of  the  roots  of  (2.  86). 

From  (2.  85),  it  is  seen  that  the  only  possibility  is  for  p2  itself  to  be  constant  and  equal  to 
ic2  .  It  can  then  be  shown  that  either  0  £,  p  =  x  <  1,  orp=  -  x  <  -  1  .  The  first  case  only  is  admissible, 
since  the  second  precludes  p  =  x  initially.  Thus,  we  have 


0ip2=  K2  -  x  i 


and 


6  >3 
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The  solution  is  than  as  follows 


1  .  k2  i  .  -2 

R  «  *  .  J_JL  vt  .  a  -  V-£-  Z 

2k  2k 


— (?) 


R\  4(*2)/(l  -  kz) 


U  =  lev 

6  =  4—  r  J5—1.  ,-t.i  >  3 

v2  2*2 


(a) 

(b) 

(c) 
(d> 


(2-87) 


For  k  =  1,  6  =  3  and  v  =  (2/'/T)c,  and  the  deformed  diaphragm  is  cylindrical,  i.e.,  it  ruptures  com¬ 
pletely. 

For  other  allowed  values  of  *  (=p),  this  solution  is  quite  similar  to  the  elementary  approx¬ 
imation  presented  earlier,  even  quantitatively,  provided  6  is  large.  For  it  follows  then 

k  « l/6l/2  =  v/2c«l,  (6»  3) 
and  (2.  87)  reduces  to  (2.  83). 


Another  case  of  interest  arises  if  k  has  any  desired  positive  value  within  reason  and  6  is 
very  large.  Then  xj  is  small  and  x2  is  large.  Assume  that  kzIx%  «  1,  as  is  \xz/x2  (since  p2  varies 
from  *2  to  xi).  Under  these  simplifying  conditions,  it  is  possible  to  derive  an  approximate  integral 
of  (2. 85  c).  It  is  then  found  that 


xj  «  1/6,  x2  «»  26  »  1 


Then  providing  6  is  large  enough 


-  =  —  =  2Uv/(v2  -  U2)  [  1  -  (R/a)] 
a  a 


(a) 

(b) 

(c) 


(2.88) 


It  may  be  shown  that  this  approximation  is  uniform  over  the  range  of  U/v  from  k  to  v/2c,  and  the 
range  of  R/a  from  1  to  0;  but  the  nearer  k  is  to  1  the  larger  6  must  be.  If  k  =  0,  which  is  the  case  of 
a  material  which  has  strictly  no  elastic  strain  range,  the  deflected  shape  of  the  diaphragm  described 
by  (2.  88  a)  is  conical  near  the  center  with  a  central  deflection  the  same  as  that  given  by  the  elemen¬ 
tary  theory  under  2.  2.  2.  2. 


2.  2.  2. 4.  General  Remarks 


It  may  be  noted  that  the  modulus  of  elasticity  E  enters  into  the  definition  of  k  .  This  is 
because,  although  the  elastic  phase  of  the  material  is  neglected  in  the  specification  of  the  deforma¬ 
tions,  an  estimate  has  been  made  of  the  initial  conditions  at  the  end  of  the  elastic  range.  The  time 
at  which  elastic  phase  is  completed  is  of  the  order  of  magnitude  of  3%  of  t(  for  some  steels.  Thus, 
it  is  generally  immaterial  whether  time  is  reckoned  from  the  end  of  the  elastic  strain  range  or 
before  it. 


Whether  the  above  analysis,  developed  with  specific  reference  to  diaphragms,  would  apply 
to  plates  of  the  proportions  of  interest  in  the  present  research  program  can  only  be  determined  by 
careful  experimentation  and  photographic  records  of  the  deformation. 


Some  additional  reaulta  which  may  be  of  value  in  checking  the  theory  are 


crr  =  irfl  =  <r 

€r  =  ee  =  log  (r/rQ) 

eH  =  log  H/h 


(2.89) 
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3.  NUMERICAL  ANALYSIS  -  CYLINDRICAL  SHELLS  (LINEAR-ELASTIC,  SMALL  DEFLECTIONS 


Th«  computation*  utilised  the  following  geometric*  end  materiel*  for  the  cylindrical  model* 
u*ed  in  the  experimental  investigation 


a  =  6" 


v  =  1/3 


L  =  36"  E  =  2.  9  X  107  (p*i) 

h  =  .036"  p  =  7.36X10*4 


The  analytes  began  with  determining  the  free  vibration  circular  frequency  utilising  the  value* 
obtained  for  the  frequency  factor  (4>)  in  Equation  (I.  13)  for  the  appropriate  values  of  k  [  see  Eqs. 

(2.  9),  (2.  16)  or  (2.  17)]  in  the  relationship 

“mn  =  <12-3W1/2xi°4  <»•« 

The  variation  of  the  circular  frequency  with  the  end-closure  stiffness  (R)  for  various  values  of  m  and 
n  are  shown  in  Figure  3.  20. 

The  computations  in  all  cases  were  directed  toward  obtaining  pertinent  radial  displace¬ 
ments  from  the  relation 


w(x.<M)  =  £wmn(x.<t>)qmn(t) 


Values  for  (t)  in  Equation  (3.  1)  were  acquired  from  Equation  (I.  27)  [  or  in  finite  difference  form, 
Eq.  (I.  28))  .’wfiich  entailed  obtaining  the  generalised  force  and  generalised  mass. 


3.  1  Simply  Supported  End-Closures 

For  the  generalised  force  [  see  Eq.  (I.  27);  Appendix  I] ,  taking  the  loading  as  being  uni¬ 
formly  distributed  along  the  length  of  the  shell,  we  have 


Q_n(t)  =  a 

mn 


f  £p(4>.  t)  f  wlnn(x,+)dxj  d* 


From  Equation  (2.  15) 


C>  (t)  =  - 

mn'  '  nit 


-  J  P($,t) 


cos  irfod$  (n  =  1 ,  3,  5,...) 


Taking  into  consideration  the  loading  symmetry  and  using  cylinder  surface  increments  of  ♦  *  22.  5*, 
Equation  (3.  4)  becomes 


Qmn^  _  n* 


4>j  +  22.  5 

P„(t)  cos  m*fo  +  X  /  P4,  +  11.  25  (t)  co*  m*d* 

♦1  1 


•  /‘ 
168.  75 


Pl80(t)  cos  m+d+ 


J 


The  value*  of  P,  (t)  for  time  increment*  of  0.  1  msec  were  obtained  from  load  profile*  aimilar  to 
those  shown  in  Figure  I.  10,  Appendix  1.  The  time  variation  of  Qj^U),  a*  acquired  from  Equation 
(3.  5)  for  Pj  *  33  pai  and  t+  *  1  msec,  is  shown  in  Figure  3. 1. 

From  Equation  (2. 15)  [  and  Eqs.  (2.  11)1 ,  the  expression  for  the  generalised  mass  [  see 
Eq.  (1.27),  Appendix  I]  becomes 

[..  (^)2*  (5£)‘] 

The  frequencies  and  the  corresponding  values  of  the  terms  within  the  brackets  of  Equation  (3.  6)  are 
given  in  Table  3.  1. 


In  each  of  the  computations,  only  the  lowest  of  three,  real  frequencies  obtained  from 
Equation  (1. 13)  was  used.  This  frequency  was  associated  with  the  response  where  the  radial  dis¬ 
placement  (w)  was  predominant.  An  indication  of  the  contributions  made  by  the  longitudinal  (u)  and 
tangential  (v)  displacements  may  be  found  in  Table  3.  1  where  sample  values  of  K^/Kj  (reflecting  the 
tangential  displacements)  and  Kjns/KjL  (reflecting  the  longitudinal  displacements)  are  given.  Here, 
it  is  seen  that,  although  these  values  are  greater  than  unity,  they  are  at  the  same  time  associated 
with  much  higher  frequencies.  Since  the  value  of  the  generalised  mass  is  increased  [  Eq.  (3.6))  and 
the  larger  frequencies  dictate  smaller  values  of  At  in  Equation  (I.  28),  the  net  result  is  a  drastic 
decrease  in  the  values  of  qmn(t).  This  clearly  shows  that  the  contribution  of  the  longitudinal  and 
tangential  displacements  is  quite  negligible  when  compared  with  that  of  the  radial  displacement. 


The  ratio  of  the  generalised  force  to  the  generalized  mass  appearing  in  Equations  (I.  27) 
and  (L  28)  is,  in  accordance  with  Equations  (3.  5)  and  (3.  6) 


0^,1.!  4  X ,/  V> «- 


(m  =  2,3,...) 
(n  =  1,3,5,...) 


(3.  7) 


where,  for  the  particular  shell  properties  previously  listed 

4  =  1.53  X  IQ6 

nhpir^  n 

From  Equation  (I.  28),  the  variations  of  qmn(t)  with  time  are  obtained  for  particular  values  of  Pj  and 
t+.  These  are  shown  in  Figures  3.  2  through  3.  13.  As  a  final  step,  Equation  (3.  2)  is  utilised  to 
obtain  the  value  of  the  radial  displacement  for  each  point  on  the  cylinder's  surface  as  defined  by  x 
and  In  Figures  3.  2  through  3.  13  are  shown  the  variation  of  w(x,  6,  t)  with  time  for  the  point  at 
x  =  L/2  and  6  =  O’.  (See  also  Table  3.  2) 

3.  2  Fixed  End-Closures 


For  cylinders  with  fixed  boundary  conditions,  the  expression  for  the  generalised  force 
[  see  Eqs.  (2.  13)  and  (3.  3)1  for  a  load  uniformly  distributed  along  the  length  of  the  shell  is 


^mnW 


1^7“  P(«)  /  P(+. »)  cos 

r  o 


where 

P(k)  =  r2(cosh  it  +  cos  k  -  2)  +  r3  ( sin  k  -  sinh  k) 
T'J  =  sinh  k  -  cosh  k  +  sin  k  +  cos  k 


(3.8) 
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FIGURE  3.  1.  VARIATION  OF  LOADING  FUNCTION  WITH  TIME 


TABLE  3.  1.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR  CYLINDRICAL 
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FIGURE  3  2.  VARIATION  OF  <Umi<*>  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  =  L/2.  4  =  0“  FOR  CYLINDRICAL  SHELL 


FIGURE  3.  3.  VARIATION  OF  qmn<t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  *  L/2,  4*0“  FOR  CYLINDRICAL  SHELL 
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#  (x,  f,  i)(mcAv  tj 


FIGURE  1.  4.  VARIATION  OF  q^lt)  AND  RADIAL.  DISPLACEMENT  w(t> 
ATi  =  L/2,  ♦  »  O'  FOR  CYLINDRICAL  SHELL 


FIGURE  1.  S.  VARIATION  OF  <L_.(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  •  L/2,  ♦  •  0*  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  6  VARIATION  OF  qmn(t)  AND  RADIAL  DISPLACEMENT  w(i) 
AT  X  :s  L/2,  <t>  -  0  FOR  CYLINDRICAL  SHELL 


T/ME  t'/wjcisJ 

FIGURE  3.  7  VARIATION  OF  qmn<0  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  L/2,  f  0*  FOR  CYLINDRICAL  SHELL 
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T/MBl  <-) 

FIGURE  3  8  VARIATION  OF  qmn(t)  AND  RADIAL  DISPLACEMENT  w(») 
AT  x  L/2.  9  =  O'  FOR  CYLINDRICAL  SHELL 


Tit*  a  (mrme*) 

FIGURE  3.  9  VARIATION  OF  <fann<t)  AND  RADIAL  DISPLACEMENT  w<») 
AT  x  ■  L/2,  ♦  •  0*  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  10:  VARIATION  OF  qWH(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  =  L/2,  ♦  «  0*  FOR  CYLINDRICAL  SHELL 


FIGURE  3.  1 1 .  VARIATION  OF  AND  RADIAL  DISPLACEMENT  w<t) 

AT  x  ■  L/2,  +  •  O’  FOR  CYLINDRICAL  SHELL 
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FIGURE  3. 12.  VARIATION  O F  w^t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  «  L/  2,  ♦  «  0*  TOR  CYLINDRICAL  SHELL 


r/MC  ) 

FIGURE  3. 13.  VARIATION  OF  q^t)  AND  RADIAL  DISPLACEMENT  w<t) 
AT  x  »  L/2,  ♦  ■  0*  FOR  CYLINDRICAL  SHELL 
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TABLE  3.  2.  COMPUTED  MAXIMUM  RADIAL  DISPLACEMENTS  (COLUMNS  A)  AND 
CORRESPONDING  TIME  OF  OCCURENCE  (COLUMNS  B)ATx>  L/2.  4>  =  0* 
FOR  CYLINDERS  WITH  SIMPLY -SUPPORTED  END-CLOSURES 


V 

Columns  A 

PI 

(1) 

(2) 

iESil 

(msec) 

(inches) 

(inches) 

17 

0.  75 

.042 

.035 

17 

2.  5 

.  150 

.  085 

23 

1.  0 

.  078 

.  065 

25 

2.0 

.  114 

.  092 

25 

3.  0 

.  152 

.  124 

33 

1.  0 

.  084 

.067 

35 

2.0 

.  164 

.  128 

37 

1.  5 

.  132 

.  107 

37 

3.  0 

.  212 

.  174 

60 

1.  0 

.  167 

.  133 

60 

2.  0 

.  284 

.  184 

60 

3.0 

.  325 

.  26 

(1)  Neglecting  contributions  of  u  &  v  in  Mmn 

(2)  Including  contributions  of  u  &  v  in  Mmn 

(3)  PlXt+/Z 


Columns  B 


(1) 

(2) 

I  (3) 

(msec) 

(msec) 

lb  -in  -msec 

0.9 

0.85 

6.  3 

1.4 

1.4 

21.  2 

0.  9 

0.9 

11.  5 

1.  3 

1.  3 

25 

1.4 

1.4 

37.  5 

1.0 

1.  0 

16.  5 

1.  4 

1.  3 

35.  0 

1.  3 

1.  3 

27.  7 

1.  3 

1.  3 

55.  5 

1.  0 

1.  0 

30 

1.4 

1.4 

60 

1.  3 

1.  2 

90 

r  2  *  sin  k  +  sinh  k 
r  j  1  co«h  K  -  COS  K 
k  *  n'e 

n'  ~  l/2(2n  +  1)  (n  =  1. 2. 3, . . . ) 


It  is  of  interest  to  note  that  the  integral  term  in  Equation  (3. 8)  is  identical  to  that  appearing  in  Equa¬ 
tion  (3.  4).  Thus,  the  values  of  the  generalised  force  for  the  simple  and  fixed  support  conditions  are 
proportional.  The  values  of  the  terms  in  Equation  (3.  8)  are  given  in  Table  3.  3. 

For  the  generalised  mass,  in  a  first  order  approximation,  the  longitudinal  and  tangential 
components  of  the  displacement  were  neglected  [  see  Eq.  (I.  27),  Appendix  il ,  so  that 


or 


,2w 


^n  “  P*  S  J  l»mn<,'*l|Z  ***♦ 


LhpawC2 
(V‘)  2k 


<4(k) 


(3.  9) 


where 


1j(k)  =  ^  —  sinh  k  cosh  k  +  k  -  cosh  x  sin  kJ 

-  cosh  2k  -  2  sinh  k  sin  k  -  j  cos  2k  J 

2/1  \ 

+  r2  I  j  sinh  k  cosh  K  -  cosh  K  sin  K  I 


The  appropriate  values  of  Oj(k)  are  also  given  in  Table  3.  3 

From  Equations  (3.  8)  and  (3.  9),  we  acquire  the  ratio  necessary  for  the  solution  of  Equa¬ 
tion  (I.  28) 


gmn<l> 

**mn 


_JL  Jkl 

phirC  oj(k) 


P^(t)  cos  m*d+ 


(m  =  2,3,...) 
(n  =  1,3,5,...) 


(3.  10) 


where  the  summation  is  identical  to  the  terms  within  the  brackets  in  Equation  (3.  5).  The  variation 
of  qmn(t)  with  time,  obtained  by  utilising  Equation  (3.  10)  in  Equation  (I.  28)  are  shown  in  Figures 
3.  14  through  3.  19.  The  variation  of  the  radial  displacement  with  time  (also  shown  in  Figs.  3.  14 
through  3.  19)  are  obtained  from  Equation  (3.  2),  or,  for  the  fixed  boundary  conditions 


„  /  Itx 

KX\ 

/  .  .  KX 

KX\ 

r  3  \c0*  TT  ' 

c°Mh  -J7  )  +  r2 

^stnh  —  - 

•in-) 

COS 

Table  3.  4  contains  the  maximum  radial  displacement  and  corresponding  times. 


12) 


In  order  to  determine  the  contributions  of  the  longitudinal  (u)  and  tangential  (v)  components 
of  displacements,  we  use  the  expression  for  generalised  mass  appearing  in  Equation  (1.27)  and  from 
Equations  (2.  11)  and  (2.  12)  obtain 
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TABLE  3.  3.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR 
CYLINDRICAL  SHELLS  WITH  FIXED  END-CLOSURES 


I 


! 

i 


m 

n 

n' 

“mn 

(rad/sec) 

f 

mn 

cycle/ aec 

PU)(1) 

I 

1 

1.5  , 

50,000 

7,960 

220 

1 

2 

2.5<3> 

1.0 

1 

3 

3.  5 

56,310 

8,960 

120,000 

2 

1 

1-5,-, 

4,640 

740 

220 

2 

2 

2.5<3) 

13,660 

2,  180 

1.0 

2 

3 

3.  5 

79, 730 

12,700 

120,000 

3 

1 

!•  5 

2,040 

325 

220 

3 

2 

2.  5(3) 

6,010 

960 

1.0 

3 

3 

3.  5 

12,050 

1,920 

120,000 

(1)  See  Eq.  3.8 

(2)  See  Eq.  3.  9 

(3)  Anti  symmetrical 

(4)  See  Eq.  3.  13 

r  (Kz\z 

W 

m 

<5)  + 

J  +  azM  \  KjL/ 

ajU)*2* 

o2(k)(4) 

k2<4> 

«1 

k3K<4> 

Kj  L 

<5) 

25,600 

11,  170 

-.  9684 

.4525 

51895 

25  X  106 

0 

-.8266 

.6167 

9.  7  X  10U 

264.8  X  1011 

-.  6338 

.6172 

25,  600 

25  X  106 

9.  7  X  10U 

11,  170 

0 

264.8  X  10U 

.  5022 

.2007 

32506 

25.  600 

25  X  106 

9.  7  X  1011 

11,  170 

0 

264.8  X  10U 

.  3335 

.0874 

28530 

45 


FIGURE  3.  14.  VARIATION  OF  ci  (t)  AND  RADIAL  DISPLACEMENT  w<t) 
AT  x  =  L II,  4  =  0*  F8k  CYLINDRICAL  SHELL 


FIGURE  3.  15.  VARIATION  OF  q(t)  AND  RADIAL  DISPLACEMENT  w(t) 
AT  x  *  L/2,  4=0*  FOR  CYLINDRICAL  SHELL 
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r/Ktm  ( 

nOUM  ).  IT.  VARIATION  OF  It)  AND  RADIAL  DISPLACE  MINT  w(t) 

All  ■  L/Z,  *  -  O'  FOR  CYLINDRICAL  SHELL 
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H  (x,  it,  s. 


(4) «/ o* 


TABLE  3.  4.  COMPUTED  MAXIMUM  RADIAL  DISPLACEMENTS  (COLUMNS  A)  AND 
CORRESPONDING  TIME  OF  OCCURENCE  (COLUMNS  B)  AT  x  =  L/2,  *  =  0* 
FOR  CYLINDERS  WITH  FIXED  END- CLOSURES 


Column*  A 

Column*  B 

PI 

*+ 

(1) 

(2) 

(1) 

(2) 

I  (3) 

(p»l) 

(m*ec) 

(inch**) 

(inches) 

(maec) 

(msec) 

(lb  -in  -m*ec) 

17 

0.75 

.030 

.025 

0.  49 

0.47 

6.3 

17 

2.  5 

.039 

.035 

0.  54 

0.  54 

21.2 

33 

1.0 

.051 

.040 

0.46 

0.  46 

16.  5 

35 

2.0 

.068 

.057 

0.60 

0.53 

17.5 

37 

1.  5 

.  049 

.  040 

0.43 

0.  48 

27.  7 

60 

2.0 

.093 

.  076 

0.  53 

0.53 

60 

(1)  Neglecting  contribution*  of  u  and  v  in  Mmn 

(2)  Including  contribution*  of  u  and  v  in 

(3) 


2 


{-.«  [*‘(^)‘]*— (^)‘J 


(3.  13) 


where 


a2(ic)  =  r  3  (  ^  sinh  *  coeh  k  +  cosh  k  sin  it)-  y  r2  Tj  (cosh  2k  +  cot  2k) 

+  r  2  ^  sinh  k  coeh  k  +  k  -  cosh  k  sin  k  ) 

Accordingly,  Equation  (3.  10)  is  replaced  by 

r„  P(k)  ^  /  P$U)  cos  m4>d«t> 


QmnW 

M _ 


phirC 


‘l<lt>  1  +  (  k7  )  +  02<k>(k^l) 


(3.  14) 


The  corresponding  values  lor  qmn(t)  obtained  from  Equation  (I.  28)  and  w(t)  at  x  =  L/2  and  4  =  0*  are 
given  in  Figures  3.  14  through  3.  19. 

3.  3  Elastic  End-Closures 

For  cylinders  with  flexible,  flat-plate,  end-closures,  we  take  as  a  first  order  approxima¬ 
tion  the  following  relation  for  the  moment 


I  s*l+ 


(3.  15) 


D  dr 

where  the  bar  denotes  properties  associated  with  the  plate.  The  solution  of  Equation  (3.  15)  for  4  is 

♦  =  “  — ~ 

D  v+  1 

Setting  the  value  of  $  at  r  =  I  equal  to  one  radian,  the  value  of  stiffness  S  is  found  by 

S  =  (for  ^  I  =1  radian) 

>  la 

If  the  flat  plate  has  the  same  thickness  as  the  cylinder  wall  and  is  made  of  the  same  material  as  the 
shell,  then  using  the  specific  dimensional  and  materials  values  previously  listed,  we  have  for  this 
particular  case 

S  =  28.  2 


and  thus 


RJ£«8 


For  a  load  uniformly  distributed  along  the  length  of  the  shell,  the  expression  for  the  gener¬ 
alised  mass  is 
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1 


Qmn(t)  =  ~  ««.*)  /  *•<♦.*»  co*  m*d* 

*  1*  ft 


(3. 16) 


where 


Ws,  R)  *  2  ^a2  [  ein  a(coah  a  -  1)  +  ainh  a(l  -  coe  a)  +  Ra(cosh  a  coa  a  -  1) 
+  R2[  ain  a(coah  a  -  1)  -  ainh  a{l  •  coa  a)]} 

=  2a2  ain  a  +  R(a  +  R)(alnh  a  -  coah  a  +  ain  a  +  coa  a) 


=  2a2  ain  a  +  R(a  +  Rjry 
r  g  =  ain  a  +  ainh  a 


(See  Eq.  (3.8)] 


r  3  =  coah  a  -  coa  a 

and  a  are  the  valuea  commensurate  with  the  value  of  R  obtained  from  Figure  2.  1.  Theae  values  of  a 
and  the  circular  frequencies  obtained  from  Figure  3.  20  are  given  in  Table  3.  5. 

For  the  generalised  mass,  neglecting  the  contributiona  of  the  longitudinal  and  transverse 
components  of  displacements,  we  have 


„  LhpawC2  , 
“nn=  - *-7—  al<*'R> 


(r-j)2  a 


(3.17) 


where 

Oj(a,  R)  =  2a®r2(ainh  a  -  ain  a)  +  2a4l  3  ainh  a  ain  a(coah  a  ain  a  -  ainh  a  coa  a) 

-  RTj  (ainh  a  -  ain  a)i  +  a3Rl  4  ainh  a  ain  a(cosh  a  cos  a  -  1)  -  Rr|] 

+  e2R2I(cosh  a  cos  a  -  l)(sinh  a  cos  a  -  coah  a  ain  a)+  8  ainh  a  ain  a  (ainh  a  cos  a 

+  coah  a  ain  a)  +  2r2rj  1  +  aR3]  RT2  +  4  (coah  a  coa  a  -  l)2) 

+  3Ri(coah  a  cos  a  -  l)(cosh  a  ain  a  +  ainh  a  coa  a) 

From  Equations  (3.  16)  and  (3.  17),  the  necessary  ratio  for  the  solution  of  Equation  (I.  28)  is  as 
follows 


Qmn<l)  rl  P(a,R)  Vf... 

^  L  I  p*«  c°' 


(3.  18) 


where  the  summation  ia  identical  to  the  terms  within  the  brackets  in  Equation  (3.  5).  The  variation 
of  w(x,4,  t)  at  x  =  L/2,  and  ^  =  0*  with  time  ia  shown  in  Figurea  3.  21  through  3.  23. 
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TABLE  3.  5.  FREQUENCIES  AND  RELATED  QUANTITIES  FOR 
CYLINDERS  WITH  EiASTIC  END-CLOSURES 


m 

2 

3 


n 

K 

u>_,_ 

mn 

(rad/eec) 

P0c.R)(1) 

OjOc.R)*2) 

1 

I.  41ir 

4,000 

12,200 

56.  49  X  1  o' 

1 

1.  41ir 

2,  100 

12,  200 

56.  49  X  1  o' 

(1)  See  Equation  {3.  16)  (for  R  =  8) 

(2)  See  Equation  (3.  17)(for  R  =  8) 
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FIGURE  3. 21.  VARIATION  Or  AND  RADIAL  DISPLACEMENT  w(t) 

AT  x  -  L/2,  t  »  0*  FOR  CYLINDRICAL  SHELL 
ii 
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FIGURE  3  22.  VARIATION  OF  qmn(tl  AND  RADIAL  DISPLACEMENT  ..-(tl 
AT  x  *  L/2.  4*0*  FOR  CYLINDRICAL  SHELL 
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FIGURE  3.  23  VARIATION  OF  qmn(l>  AND  RADIAL  DISPLACEMENT  «•('! 
AT  x  -  L/2,  (>0>  FOR  CYLINDRICAL  SHELL 


4.  DISCUSSION 


The  utility  of  an  analytical  procedure  ia  beet  meatured  in  terma  of  ite  ability  to  yield 
meaningful  and  correlative  quantitiee  for  a  wide  range  of  variablee  with  a  minimum  of  computational 
effort.  The  adequacy  of  theee  computed  quantitiee  can  be  aecertained  only  in  a  direct  comparieon 
with  reliable,  equivalent  experimental  evidence.  If  the  reeulta  obtained  by  a  comparatively  aimple 
analytical  method  agree  favorably  with  the  available  experimental  data  and  if  the  reeulta  cover  the 
range  of  intereet  for  the  problem  at  hand,  there  then  appeara  to  be  little  reaeon  (at  leaat  from  a 
practical  viewpoint)  for  employing  another  more  rigoroua  and  therefore  more  computationally  com¬ 
plex  eolution. 

The  linear- elaetic  method  for  cylindrical  ahella  developed  in  Section  2.  1. 1  and  used  in  the 
numerical  analyeis  in  Section  3  appeare  to  yield  eignificant  results  (over  a  wide  range  of  peak  over¬ 
pressures  and  positive  phase  durations)  that  agree  favorably  with  the  experimental  data.  There  is, 
at  present,  no  apparent  or  justifiable  reason  for  suggesting  the  use  of  other  available  but  neverthe¬ 
less  more  complex  analytical  procedures  (such  as  that  developed  in  Section  2. 1. 2). 

One  method  of  compactly  presenting  the  numerical  results  obtained  in  the  preceding  sec¬ 
tion  is  to  chart  the  variation  of  the  maximum  radial  displacement  (in  this  case,  at  x  -  L/2,  and 
4=0*)  with  the  two  characteristics  which  distinguish  each  particular  blast  loading,  namely,  the 
peak  incident  overpressure  (Pj)  and  the  companion  positive  phase  duration  (t+).  In  Figure  4.  1,  we 
have  taken  the  impulse  (defined  as  Pj  X  t+  -r  2)  for  the  ordinate  axis  and  the  maximum  radial  deflec¬ 
tion  for  the  abscissa.  On  one  hand,  the  bands  for  each  type  of  cylinder  end-closure  arise  from 
taking  the  generalized  mass  as  a  function  of  only  the  radial  component  of  displacement  and  on  the 
other  including  all  three  displacement  components  (u,  v  and  w)  in  computing  Mmn.  In  Table  4.  1  are 
comparisons  of  high-speed  camera  data  obtained  for  specific  models  with  the  computational  results 
as  represented  in  Tables  3.  2  and  3.  4  or  Figure  4.  1. 

The  close  agreement  between  the  experimental  and  analytical  results  is  indeed  encour¬ 
aging  particularly  in  view  of  some  of  the  expected  deviations  in  the  peak  incident  overpressure  and 
positive  phase  duration  that  arise  from  minor  differences  in  charge  size  and  shape,  differences  in 
atmospheric  conditions  and  local  changes  at  ground  zero  resulting  from  repeated  shots.  Moreover, 
there  are  the  assumptions  used  in  developing  the  load  spectrums  and  the  purely  analytical  approxi¬ 
mations  (see  Appendix  I)  and  some  of  the  normal  uncertainties  associated  with  the  high  speed  photog¬ 
raphy  data  (see  Appendix  III)  and  the  models'  dimensional  and  materials  properties.  All  these  fac¬ 
tors  undoubtedly  contribute  to  the  development  of  a  range  of  displacement  values  rather  than  a 
precise  and  specific  value  for  each  given  set  of  conditions.  The  encouraging  aspect  is  that  the 
numerical  procedure  used  in  the  preceding  section,  which  is  comparatively  simple  and  straight¬ 
forward,  does  give  values  that  are  compatible  with  the  experimental  evidence  and,  therefore,  of 
direct  practical  importance. 
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FIGURE  4.  1.  VARIATION  OF  MAXIMUM  RADIAL  DISPLACEMENT  WITH  ONE 
HALF  OF  PRODUCT  OF  PEAK  INCIDENT  OVERPRESSURE  AND 
POSITIVE  PHASE  DURATION 
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TABLE  4. 1 .  COMPARISON  OF  EXPERIMENTAL  (HIGH-SPEED  PHOTOGRAPHY) 
AND  COMPUTED  (LINEAR- ELASTIC,  SMALL  DEFLECTION  THEORY) 
RADIAL  DISPLACEMENTS  FOR  CYLINDRICAL' SHELLS 


Model 

No. 

Shot 

No. 

Boundary 

Restraint 

*1 

(P»i) 

*+ 

(msec) 

(PlXtt)/2 
(lb  -  in1 2 3 4 5 6 7 8 9 10  -  msec) 

Computed 
(5)  ((>) 

(inches)  (msec) 

Experimental'^' 

(Ti  w 

(inches)  (msec) 

CA3 

D1 

Fixed 

15(2) 

6.  3  (4) 

47 

.064  -  .078  (7) 

. 

.04 

0.7 

CA3 

D2 

It 

19  (2) 

5.  5  (4) 

41 

.060  -  .072  (7) 

- 

.07 

1.  5 

CA3 

D3 

II 

42  (2) 

3.8(4) 

80 

.098  -  .  118  (7) 

\ 

.  12 

0.8 

CB2 

G1 

8.  8. 

20  (1) 

2.4(3) 

24 

.084  -  .  116  (7) 

1.4(10) 

.  12 

1.3 

CB2 

G2 

it 

16  (1) 

2.3  (3) 

19 

.072  -  .095  (7) 

1.4(10) 

.  11 

1.  5 

CB2 

G3 

tl 

17  (1) 

2.5(3) 

21.5 

.085  -  .  105  (8) 

1.4(8) 

.  11 

1.4 

CB2 

G4 

it 

17  (2) 

2.9(3) 

24.  5 

.088  -  .118  (7) 

1.4(10) 

.  13 

1.3 

CB2 

HI 

H 

25  (2) 

5.  1  (3) 

64 

.194  -  .248  (7) 

- 

.  20 

2.0 

CB3 

H3 

It 

23  (1) 

4.5(3) 

52 

.164  -  .216(7) 

- 

.  15 

1.5 

CB3 

H4 

II 

53  (1) 

3.9  (3) 

100 

.  284  -  .  366  (8) 

1.4(8) 

.  30 

2.0 

(1)  From  Table  III.  1,  Appendix  III 

(2)  From  Figure  in.  3,  Appendix  III 

(3)  From  Table  III.  2,  Appendix  III 

(4)  Computed 

(5)  Maximum  (inward)  radial  displacement  at  x  =  L/2,  4>  =  0' 

(6)  Time  of  occurrence  of  displacement 

(7)  From  Figure  4.  1 

(8)  From  Table  3.  2  or  computed 

(9)  From  Figures  III.  4  through  III.  13,  Appendix  III 

(10)  Estimates 


APPENDIX  I 

METHODS  OF  ANALYSIS  (UNEAR-ELASTIC. 
SMALL  DEFLECTION  THEORY) 


LIST  OF  SYMBOLS 


A.  B,  C 


Constant  coefficient*  (maximum  amplitudes  of  component  vibrations) 


12(1  -  v  ) 
Elastic  modulus 
pa2(l  -  v2) 


Cylinder  length 
Bending  moment  of!  section 


Generalised  mass 


Stress  resultant 


Generalised  force 


Kinetic  energy 
Potential  ( strain)  energy 
Mean  cylinder  radius 
Cylinder  wall  thickness 


Subscript 


Number  of  circumferential  waves 


Generalised  coordinate 


Longitudinal,  tangential  and  radial  components  of  displacements 
Axial,  tangential,  and  radial  coordinates 


I 


I 


LIST  OF  SYMBOLS  (Cont'd) 


a  =  ^/L 

y  Shear  etrain 

t  Axial  (train 

v  Poisaon's  ratio 

(lb- sec*  \ 
— ~ —  / 

1  -  v2  22 

4,  =  — g —  pu  a 

oj  Circular  frequency 


I 


V 

l 
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1.1  FREQUENCY  AND  DISPLACEMENT  EQUATIONS 


The  development  of  the  necessary  free-vibration  equations  follows  the  format  estab¬ 
lished  by  Yu.  (I*  1)  The  solution  begins  with  the  following  equations  of  motion(  !•  2)  based  on  the 
assumption  that  N^r  =  Nxr  *  0.  (See  Figs.  1. 1  and  I.  2  for  notation.) 


8N  8N 


a  -  pahU  =  0 


»N  .  8N 

a  TT  +  ST  +  pahV  =  0 


»2“*+  »2mxx  1  »2M»x 


^+T^?-+-?CTr+i-^5-  +  N<H-pabW  =  0 


(I.D 


Based  on  the  condition  that  a»z  (so  that  a  +  z  -  a)  and  the  following  equatione  for  strains  as 
functions  of  displacements 


8u  8^w 


SH 


8x2 


1  3v  w  z  9  4w 
'  I  84-  +  a  '  “2  172 " 


3v  .  1_  8v_  s  8  w 
a  8+  " C  a  8x8<|) 


(1.2) 


we  obtain  the  following  differential  equations  for  the  displacements 


4  1  83w 

V  u  +  I  - 2 

a  8x8y^ 


v  83w  F  D2  |3.»  2  2v  8w 

*a^-7^L~  117775  57 


2F  8 


a  (1  -  v)  8t 


>]■ 


4  2 

V4v  -  - 


+  V  83w  _1_  8^W  F  82  [ 

*  8x28y  *  8y3  a2  8t2  [ 

2F  82v  1 

(1  -  v)  it2’] 


3  -  v  v2y 


a(l  -  v)  8y  .2 


(1.3) 


d.4) 
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(1.5) 


D  _8  1  -  v2  84w 

K  V  w  + 5 - T 

K  a2  8x4 


_F  82 
a*  dt^ 


[(- 


2F 


a2(l 


8\  /  F  8  i 

v)  8t2'  '7  7? 


w 

*7* 


-  -* — 
*2(1 


7,( 


82w 


2  8Zw 
w  — t  +  — r 

8x2  8y2 


)] 


where 


y 

F 

D 


=  a<> 

pa2(l  -  v2) 
E 

Eh3 

'  12(1  -  v2) 


K 


The  space  and  time  dependent  forms  of  the  displacement  equations  are  taken  to  be 


where 


UjnnU.+.t)  =  £  Ajexaj  cos  m*  sin  <*nnt 

j 

vmn(x>  *)  ■  Z  BjeXai  »»n  m$  sin  u*nnt 

j 

wmn(x'  *1  =  Z  Cj«XOj  co»  ">+  sin  <^,,1 
j 


r  (j  =  1.2 . 8> 


«{  =  -* 


(1.6) 

(1.7) 

(1.8) 


m  *  number  of  circumferential  waves 

(2m  =  number  of  circumferential  nodes) 

Substituting  Equations  (1.  6),  (1.7),  and  (1. 8)  into  Equations  (I.  3),  (1.4),  and  (I.  5),  we  have 

?  “/“'{rrr  frrVX  [>  •  (-£-)  ]  ♦«*  [‘  -(-£)]} 

2  [■  -  (2  *  v)  (£)  ]} 


+  m 


=  0 


(1.9) 
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(1. 10) 


£  +  frr  w>2  [i  -  ®  ]  +  [*  -  (¥)]} 

-  ajCje^j  jm2  [l  +  v  Q-)  ]  -  =  0 

Z  Cj®xaj  [,L .  ^n|2  +  (3  -  v)m2  [i  . 

+  Vn{(‘  -  ‘')«2  [l  -  (l?)  ]  +  <1  -  Wm4  [l  -  [-£)  ] 

‘  [.-(5^-c..^ [..(?)] 

2)(aaj)4j  =  0 


+  (3  -  w)kmc 


where 


(1  -  v)(l  -  v‘ 


1  -  v  2  2 

l-PVn* 


12a‘ 


Equations  (I.  9).  (I.  10),  and  (1. 11)  are  simplified  by  assuming  that  the  shell  wall 
thickness  (a)  and  length  (L>)  are  such  that 

-*  -.v2 


(-0  -  Q)  «>■• 


This  simplification  permits  Equations  (I.  6),  (I.  7),  and  (I.  8)  to  be  rewritten  as  follows 


umn<x,4’-t>  *  iq-  cos  *"♦  "in  “mnl  Z  ajCjexaj  1 


K, 


vmn<x-<t>-  »)  =  k7  sin  m4>  sin  <%mt 


zCj.“j  L 


<t>.  t)  =  cos  sin  u^t  L,  Cje  J 

j 


(j=l,2,3.4) 


(Ml) 


(I-  12) 
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where 


K1  =  ^mn  *  <3  *  v>  W"2  +  <*  *  v>m4 
K2  =  m{m2(l  -  v)  -  24<mnJ 
K3  =  «(2*mnv  +  (1  -  v)  m2J 


From  Equation  (I.  11),  we  also  obtain  the  following  frequency  equation 

^mn  *  'J'mnl2  +  (3  -  v)m2  +  2km4]  +  ’('mnK1  -  w>m2  +  (1  -  v)m4 

+  (3  -  v)km6l  -  (1  -  v)km8  -  (1  -  v)(l  -  v2)(aaj)4  =  0  (j=l,2,3,4)  (1.13) 

The  solution  of  Equation  (I.  13)  for  4*mn  (and  therefore  the  frequency  Ug,,,)  depends 
on  obtaining  the  appropriate  values  of  Uj(j  =  1,2,3,  4),  which  reflect  the  shell's  end-closure 
boundary  condition.  For  the  fourth  degree  equation,  the  expected  values  of  Uj  would  be  of  the 
form 


s  K  is  is 

“1  =  L  •  a2  =  -  l  •  a3  =  T  and  a4  =  *  T 


where  i  =  «/TT  and  s  is  a  real  number  containing  a  term  reflecting  the  axial  mode  shape. 

I.  2  Lagrangian  Equations  of  Motion 

The  development  of  the  Lagrangian  equations  of  rpotion  necessary  in  the  numerical 
analyses  for  the  response  of  cylindrical  shells  to  impulsive  loading  begins  with  the  general  dis- 


placement  equations 

u(x.+,  t)  = 

^umn(x'+'t) 

v(x,  4>,  t)  = 

Z  vm„(x.  ‘) 

► 

(X.M) 

w(x.+,t)  = 

£wmn<x-+'‘> 

Utilizing  Equations  (I. 

12),  these  may  be  rewritten  in  terms  of  generalized  coordinates 

u(x,  <M)  = 

^  “mnl*1 

v(x.<M)  = 

Z  vmn(x’  tllmnl1! 

* 

(1.15) 

w(x,  <>,t)  = 

Z  wmn(x, 

The  expression  for  the  kinetic  energy  is 


=  £ 
2 


t  /'  J 


h/2 


[u(x,  4».  t)]‘ 


+  l  6<x,  $,  t))  2  + Wx,  4>,  *)]  2  ]■  adxd^dz 


1' 


(1. 16) 


Since  6(x,  t) 
be  rewritten  as 


Z  umn(x*  $)4mnW  (and  similarly  for  v  and  w),  Equation  (I.  16)  may 


T  =  X  f  £  {[umn(x.«)4mn<t)J2  +  Ivmn<x,+>Wt>l2 

+  fwmn<x*  O’mnWI  2  +  “mi>(x'  ♦Kjt**  4>)4mn(t>4ij(t) 

+  VJX'  ♦>vij<3C*  ♦> +  Wmr>1*'  ♦)Wij(x*  ♦MV*1}  ****  (L  17) 


where  mn  4  ij 


Since  the  generalized  coordinates  are  the  principal  or  normal  coordinates,  the  corre¬ 
sponding  vibrations  are  the  principal  modes  of  vibration,  and  the  products  of  the  velocities  in 
Equation  (1.  17)  vanish.  That  is,  from  Equations  (1.  12),  any  one  of  the  last  three  terms  in 
Equations  (I.  17)  would  be  of  the  form 


Wij  /  L  Fm(x)Fi(x>  /  2"  Gn(4)Gj($)d$ 

X  =  o  l<*>  =  o 


dx 


where  Gn($)  =  sin  or  cos  n$,  and  Gj(4)  =  sin  or  cos  j$. 
Thus 


*=  2ir 


/  Gn(<l>)Gj($)d<t>  =  0 


6  -  0 


(for  n  4  j) 


and  Equation  (1.  17)  becomes 


T  =  ^  Z  lv(‘»Zl$umn  +  *mn  +  "mnl 

where 

Vn  =  f  f  lumn(x,6)]2  dxd<> 


(1.  18) 


(I.  19) 


and  similarly  for  vmn  and  wmn. 
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The  expression  for  the  potential  (etrain)  energy  in  terms  of  the  derivative  of  the  dis¬ 
placement  components  and  the  generalised  coordinates  may  be  obtained  by  utilizing  Equations  (I.  2) 
and  (1. 15)in  the  following 


Ea  f+h/2  p  1  -  v  2  1 

- T~  J  J  J  ef  +  ef  +  2vc„Ci  +  —5— 

2d  -  O  -h/2  0  0  L  *  *  X*  2  *XJ 


dxd^dz 


(120) 


Since  terms  containing  products  of  the  generalized  coordinates  vanish,  Equation  (1.  20) 
may  be  rewritten  as  follows 


V  = 


Ea 


2(1  -  v4) 


llmn(*)l  3  [h/T  F(u,  v,  w)dxd$^ 


h3  1 

~T  J  J  G(u,  v,  w)dxd4> 

0  0  J 


(1.21) 


where  F(u,  v,  w)  and  G(u,  v,  w)  denote  a  collection  of  terma  involving  the  partial  derivatives  of 
%„(*.♦).  vmn(x,  $),  and  *mn(z,  *)  with  respect  to  x  and  +. 

From  Equation  (I.  21),  we  have 


8<wXt)  *  H<u'v'w>Wt>  (I*2> 

where  H(u,  v,  w)  are  merely  the  combined  constant  and  integral  terms. 

For  Lagrange's  equations,  we  have 

d  T  8T  1  8T  8  V  _  ... 

<«  ‘  +  ^  =  ^"n( 

From  Equation  (I.  20) 


8T 

■ 0 

d  I  8T  I  *  .*  t 

dt  I  84mn(t)  I  =  ^W^Snn  +  vmn  +  wmnl 

Introducing  Equations  (I.  22),  (1.  24),  and  (I.  25)  into  Equation  (I.  22) 
P^Snn  +  *mn  +  *mnl  +  H<u’ v' w>  W*>  = 


(I.  24) 


(125) 


(126) 
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Expressing  Equation  (I.  26)  in  a  more  convenient  form,  we  have 


W*1  +  “mnVn(t>  =  TT““ 


where 


“mn  =  P^l^mn  +  *mn  +  ^mnl  (generalised  mass) 


H(u,  v,  w) 
M _ 


Qj^ft)  =  Generalized  force 
.L  ~2ir 


■u 


P(x,  <(».  t)  w  (x,  6)  adxd<fr 


(I.  27) 


For  the  numerical  solution  of  Equation  (I.  27),  the  method  of  finite  differences  was 
used,  and,  taking  the  second  central  differences,  we  have 


'Jmn^r)  *  2<lmn^ti^  +  ^mn^f) 
_t? 


or 


Qmn<V 


(At)  +  [2  -  wmr^At)  ]  qmn(tj)  -  qmn(tjf ) 


(1.28) 


where 


tr  =  tj  +  At 

t,  ’=  t.  -  At 
I  l 


I.  3  Development  of  Load  Spectrums 

In  obtaining  the  expression  for  the  generalized  force  in  Equation  (I,  27),  use  was  made 
of  the  loading  information  for  cylindrical  shells  contained  in  References  I.  3  and  I.  4,  as  well  as 
other  publications.  Basically,  the  time  variation  of  the  loading  along  any  one  element  of  the 
cylinder  can  be  considered  of  the  form  shown  in  Figure  1.  3  or  I.  4,  the  former  utilizing  the  infor¬ 
mation  from  Reference  1.  3,  and  the  latter  from  Reference  1.  4. 

In  both  cases,  time  is  measured  from  the  instant  the  shock  front  strikes  the  leading 
edge  (element  at  6  =  0*)  of  the  cylinder.  Again,  both  forms  include  the  arrival  time  (ta)  for  each 
particular  cylinder  element,  the  duration  of  the  diffraction  loading  phase  (tjj)  and  the  positive 
phase  duration  (t+).  The  drag  phase  (for  ?  >  tD)  is  given  by 

P  =  Pj(?)  +  Cd(6)q(?)  (1.29) 
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FIGURE  I.  4.  SCHEMATIC  OF  PRESSURE  PROFILE  FOR  SHELL 
AT  ELEMENT  4-4, 


where,  for  the  time  increment  0  <  ?  <  t+  measured  chronologically  from  ta,  the  relationship  for 
Pj  (?)  and  q(t)  is  given  in  Figure"!.  57  and  C<}(+)  is  the  drag  coefficient  for  the  particular  cylinder 
element  under  consideration.  For  the  loading  shown  in  Figure  I.  3,  the  drag  coefficient  0^(4)  is 
obtained  from  Figure  6.82b,  page  272,  Reference  1.3;  for  the  loading  shown  in  .Figure  1.4,  C,j(4) 
is  found  from  the  curve  reproduced  in  Figure  I.  6. 

The  essential  difference  between  the  pressure  profiles  shown  in  Figures  I.  3  and  I.  4  lies 
in  the  diffraction  phase.  Since  the  pressure  profile  in  Figure  1.  4  represented  a  more  realistic 
loading  condition,  there  was  reason  to  tend  towards  this  approach.  Taking  into  acccount  the  experi¬ 
mental  model's  dimensional  ratios  (L/D),  the  crossover  times  (D/U)  and  the  anticipated  (as  well  as 
experimentally  recorded)  durations,  there  appeared  to  be  adequate  aerodynamic  similitude  to 
justify  using  the  experimental  curves  in  Reference  1.4  in  determining  the  appropriate  values  of 
Pj,  tj  and  tp  for  this  investigation's  loadings.  Moreover,  the  maximum  deformations  were  found 
to  occur  at  times  approximately  equal  to  or  less  than  the  diffraction  phase  duration  (tp).  This 
meant  that  the  diffraction  phase  constituted  the  major  (and,  in  some  cases,  the  entire)  deforma- 
tional  loading. 

The  curves  used  in  determining  the  significant  values  of  time  and  loading  required  to 
define  the  load  spectrum  in  Figure  I.  4  are  given  in  Figures  I.  7  through  I.  9.  The  basic,  computa¬ 
tional  procedure  was  as  follows:  the  values  of  the  peak  incident  overpressure  (Pj)  and  positive 
phase  duration  (t+)  were  obtained  from  Tables  III.  1  and  III.  2.  The  variation  of  the  incident  over¬ 
pressure  and  the  dynamic  pressure  (see  Fig.  I.  5)  were  found  from  the  following  relations 


(I.  30) 

(1-31) 

where 

2.  5Pj 

=  7P '  tp:  <P8i> 

o  X 

PQ  =  Ambient  pressure 

Next,  the  crossover  time  (D/U)  was  computed,  where 


Pl(t) 

’  pi( 

:-*) 

“p  (-  t) 

q(t) 

■  id 

M: 

)  -(•!)] 

(ft  per  sec) 


(I.  32) 


CQ=  Ambient  sound  velocity  (ft  per  sec) 


From  Figures  I.  7  and  I.  8,  the  appropriate  values  for  ta>  tj  and  tp  were  obtained  for  each  element 
in  the  cylinder's  surface.  A  numerical  analysis  to  determine  the  optimum  number  of  elements  to 
be  considered  showed  that  for  increments  of  =  10",  15*  or  22.  5*,  only  minor  differences 
existed  in  the  final  numerical  results.  Thus,  in  all  subsequent  computations,  the  elements  con¬ 
sidered  were  those  at  4  =  0*,  22.5*,  45*,  ....  180*. 

The  values  of  Rj(4)  for  each  element  were  obtained  from  Figure  I.  9.  The  appropriate 
value  of  Pj  corresponding  to  the  time  tj  was  determined  from  Equation  (I,  30)  (where  t  =  tj)  or 
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if! 

from  Pj(t)  curve*  *uch  a*  are  shown  in  Figure  I.  5.  The  critical  pressures  were  then  found  for 
each  element 


Pi  =  *i*«pi*tl> 


Typical  values  for  critical  pressures  and  times  are  given  in  Table  I.  1;  the  corresponding 
load  spectrums  are  shown  in  Figure  1.  10. 
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FIGURE  I.  6.  VARIATION  OF  DRAG  COEFFICIENT  WITH  CYLINDER  ELEMENT 
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FIGURE  I.  7.  TIME  OF  ARRIVAL  OF  SHOCK  FRONT  ON  CYLINDER  SURFACE 


DURATION  (tD)  WITH  CYLINDER  ELEMENT 
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FIGURE  L  7.  TIME  OF  ARRIVAL  OF  SHOCK  FRONT  ON  CYLINDER  SURFACE 


FIGURE  I.  8.  VARIATION  OF  CRITICAL  TIMES  (tj)  AND  DIFFRACTION  PHASE 
DURATION  (tD)  WITH  CYLINDER  ELEMENT 
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FIGURE  I.  9  VARIATION  OF  CRITICAL  PRESSURE  (Pt)  TO  PEAK  INCIDENT  OVERPRESSURE  (Pj) 

RATIOS  WITH  CYLINDER  ELEMENT 
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FIGURE  I.  10.  PRESSURE  PROFILES  ALONG  VARIOUS  ELEMENTS  OF  CYLINDRICAL  SURFACE 


TABLE  I.  1.  CRITICAL  TIMES  AND  PRESSURES  (Pj  =  33  psi 
t+  ■  1  miec)  D/U  *  0.  S3  msec 
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APPENDIX  n 


EXPERIMENTAL  PROGRAM 


DU.  OBJECTIVES 

It  was  da •  mad  advisable  to  undartaka  a  llmitad  axparimantal  program  which  would  provida 
data  to  taat  tha  adaquacy  ol  tha  analytical  procaduraa.  Knowing  that  tha  analytical  rasulta  would  ba 
in  tha  form  of  pradictad  displacements,  an  axparimantal  technique  waa  avolvad  for  obtaining  racords 
of  tha  raaponsa  of  cylindrical  ahalla  or  flat  plataa  to  known  blast  loadings . 

II.  2.  MODEL  DESIGN 

11.2.1  Cylindrical  Shalla 

All  the  cylindrical  shall  modal*  had  tha  same  overall  nominal  dimensions:  diameter  -  12", 
length  -  36",  shall  wall  thickness  -  0.  036".  The  three  baaic  configurations  differed  with  respect  to 
the  and-cloeures  and  the  related  shell  boundary  condition.  For  the  Type  CA  specimen  (Figure  11.2), 
the  shell  was  welded  to  a  1/2"  thick  end  plate  so  as  to  closely  approximate  a  shall  with  fixed-end 
conditions.  Tha  end-closures  for  the  Type  CB  specimen  simulated  a  shell  with  a  simply- supported 
boundary  condition  (Figure  II.  3).  Type  CC  shell  models  use  end  plates  of  the  same  nominal  thick¬ 
ness  a*  the  shell  (0.036").  These  plate*  were  welded  to  the  shell  as  is  shown  in  Figure  II.  1.  The 
Type  CC  models  provided  an  elastic  boundary  and  ware  used  to  demonstrate  the  effect  of  end 
restraints  of  flat  plate  end- closures. 

Within  each  CB  type  model,  a  scribe  plate  was  mounted  at  mid- length  on  the  three  support 
rods  (see  Figure  H.  3).  Four,  spring-loaded  styli,  supported  on  the  scribe  plate  and  bearing  against 
the  shell  wall,  were  used  to  record  the  deflection  at  mid-plane  points  on  the  0*,  90*,  180*,  and  270* 
elements  of  cylinder  (0*  denoting  the  blast  sida  of  the  cylinder).  A  black  longitudinal  and  circumfer¬ 
ential  grid  on  one-inch  spacing  was  painted  on  the  yellow  surfaces  of  the  models  used  in  series  D,  A, 
G  and  H. 


II.  2.  2.  Circular  Flat  Plates 


The  test  fixture  used  to  support  the  circular  flat  plates  is  shown  in  Figure  II.  4.  As  shown 
in  the  drawing,  the  plate  is  simply  supported  at  the  boundary;  for  the  fixed  boundary  condition,  an 
oversised  plate  was  used  and  clamped  between  the  inner  pipe  and  the  outer  ring.  In  order  to  evaluate 
the  effects  of  membrane  and  beam  action,  0.0312"  and  0.0598"  thick  plates  were  tested.  White 
radial  lines  at  45*  were  painted  on  blue  plate  surfaces. 

II.  3.  INSTRUMENTATION 


Instrumentation  consisted  of  synchronised,  high-speed  photography  to  record  the  cylinders' 
dynamic  response,  mechanical  styli  recordings  of  maximum  deformations,  and  a  quartz  transducer 
system  to  record  overpressures  and  positive  phase  duration*.  Each  model  was  covered  by  one  or 
two  high-speed  cameras  with  film  speeds  from  3,  000  to  5,  000  frames  per  second.  A  sequence  cir¬ 
cuit  was  constructed  to  synchronise  the  detonation  and  camera  action  so  that  the  initial  impact  of  the 
shock  wave  on  the  model  would  occur  when  the  cameras  had  obtained  their  maximum  speed. 

Two  quarts  transducers  were  flush-mounted  in  a  18"  X  24"  steel  surface  plate.  An  adap¬ 
tor,  consisting  of  a  double  brass  fitting  separated  by  a  viscous  potting  material,  was  constructed  to 
minimise  ground  acceleration  effects  in  the  pickup.  The  transducers  were  placed  at  ground  level  at 
radial  distances  from  ground  zero  equal  to  those  from  ground  zero  to  the  model.  The  transducers' 
outputs  were  fed  through  arhplifier- calibrators  to  an  oscilloscope  and  recorded  by  oscillograph 
cameras. 

II.4.  TEST  PROGRAM 


A  total  of  nineteen  charges  in  nine  eerie*  was  detonated.  The  charge*  used  and  the  cylin¬ 
drical  shell  model*  and  flat  plates  exposed  in  each  shot  are  given  in  Table  U.  1.  Each  cylindrical 
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FIGURE  II.  1.  TYPE  CC  CYLINDRICAL  SHELL  MODEL 
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TABLE  H.  1 .  SUMMARY  OF  EXPERIMENTAL  PROGRAM 
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Not  useable  None  Good 


■hall  model  was  suspended  with  nylon  cords;  the  position  and  orientation  of  each  model  relative  to 
ground  aero  are  shown  in  Figures  II.  5  through  II.  8;  the  controlling  distances  are  given  in  Table  II.  1 . 

Shot  series  A  through  C  were  used  to  obtain  qualitative  information  on  the  deformational 
response  of  the  cylindrical  shells  as  influenced  by  such  factors  as  the  distance  and  orientation  of 
models  with  respect  to  ground  zero.  In  particular,  the  tests  were  used  to  determine  the  modes  of 
failure  (in  the  form  of  permanent  deformation  or  fracture)  as  related  to  estimated  overpressure. 

The  tests  using  5  lb  of  flaked  TNT  (Series  E  and  G)  were  conducted  on  Southwest  Research 
Institute  grounds  for  the  purpose  of  checking  out  the  equipment  and  instrumentation.  In  Series  G,  the 
distance  between  the  model  and  the  charge  was  10  feet,  and,  at  this  distance,  it  was  possible  to 
photograph  small  (but  measurable)  elastic  deformations  and  mode  shapes.  The  ground  zero-to- 
model  distances  for  the  30-lb  shots  at  Camp  Bullis  varied  from  Zb  feet  to  15  feet;  each  series  of 
■hots  terminated  with  failure  of  the  model  in  the  form  of  large,  permanent  deformations.  (See 
Figures  II.  9  through  II.  1 7 . ) 
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FIGURE  II.  5.  MODEL  POSITIONS  FOR  SHOTS  Al,  A2,  AND  A3 


FIGURE  II.  6.  MODEL  POSITIONS  FOR  SHOTS  Bl,  B2,  AND  B3 
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FIGURE  II.  7.  MODEL  POSITIONS  FOR  SHOTS  Cl  AND  C2 
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FIGURE  II.  8.  MODEL  POSITIONS  FOR  SHOTS  D,  F,  G,  AND  H 
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i  FIGURE  II.  9.  END  PLATE  FAILURE  IN 

MODEL  CC1  AFTER  SHOT  A3 

i 
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FIGURE  II.  10.  END  PLATE  AND  SHELL  FAILURE 
IN  MODEL  CC1  AFTER  SHOT  B3 
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FIGURE  II.  1 1.  SLIGHT  REDUCTION  IN  SHELL  DIAMETER 
IN  MODEL  CC2  AFTER  SHOT  B1 


FIGURE  II.  12.  SHELL  AND  END  PLATES  PERMANENT 
DEFORMATIONS  IN  MODEL  CC2  AFTER  SHOT  B2 


FIGURE  II.  13.  SHELL  AND  END  PLATE  PERMANENT 
DEFORMATIONS  IN  MODEL  CC3  AFTER  SHOT  Cl 


FIGURE  II.  14.  END  PLATE  DEFORMATIONS  IN  MODEL  CC4 
AFTER  SHOT  C7 
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FIGURE  n.  15.  END  PLATE  DEFORMATION  IN  MODEL  CCS 
AFTER  SHOT  8 


.  SHELL  DEFORMATION  IN  MODEL  CA2 
AFTER  SHOT  7 


FIGURE  II.  16 
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FIGURE  U.  17.  FAILURE  STATE  OF  MODELS  USED  IN  SHOTS  D,  F,  AND 
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APPENDIX  Ill 


EXPERIMENTAL  RESULTS 
m.  1.  PRESSURE- TIME  MEASUREMENTS 


In  developing  the  experimental  portion  of  the  program,  consideration  was  given  to  creating 
blast  phenomena  with  conventional,  chemical  high  exploeives  which,  when  applied  in  an  analysis  pro¬ 
cedure  to  scaled  models  of  missile  bodies,  would  provide  a  discernible  equivalence  for  the  nuclear 
weapon  effects  on  actual  missile  structures.  To  this  end,  full  use  was  made  of  the  information  avail¬ 
able  in  Reference  III.  1,  as  well  as  other  supporting  publications  containing  information  on  HE  explo¬ 
sions. 


The  peak  shock  overpressure,  as  related  to  distance,  plays  a  leading  role  in  explosion- 
damage  correlations.  Selection  of  the  controlling  parameters  for  the  HE  blast  phenomena,  there¬ 
fore,  began  with  the  relation  shown  in  Figure  III.  1  between  the  peak  overpressure  in  free  air  and 
scaled  (radial)  range  as  proposed  by  the  Kirkwood- Brinkley  theory(EL  2)  substantiated  by  the 
experimental  evidence  of  Fisher^- 3)  and  Weibull(^-  In  this  relation,  the  scaled  range  is 
defined  as 


X  * 


R 

(W)l/3 


where 


R  is  the  radius  (feet) 

W  is  the  weight  (yield)  (pounds) 


Since  the  HE  explosions  would  take  place  on  or  near  the  ground,  other  blast  phenomena 
associated  with  surface  or  near  surface  bursts  were  taken  into  account.  Foremost  among  these 
considerations  was  the  ground  reflection  of  the  blast  wave  and  the  subsequent  coalescence  of  the 
incident  and  reflected  waves  in  the  development  of  the  Mach  stem. 


In  accordance  with  one  definition,  let  Wr  be  the  sise  of  burst  exploded  over  a  reflector  at 
a  height  of  burst  equal  to  hR.  At  a  slant  range  of  Rr  >  2  hj  (thereby  placing  the  point  in  question  in 
the  faV  Mach  region),  the  reflected  peak  overpressure  is  Pr.  For  a  free-air  burst  of  Wj,  the  inci¬ 
dent  (or  free-air)  peak  overpressure  at  a  range  of  Rj  equal  to  Rr  would  be  Pj.  Let  the  yield  of  Wj 
be  such  that,  at  a  range  of  Ri  =  Rr,  the  incident  and  reflected  peak  overpressures  are  equal.  Since 
the  scaled  slant  ranges  are 

X-R=(^73  Xi=^T73 


by  definition 

Rp  (reflection factor)  =  — i- 
WR 


3 


Accordingly,  for  a  given  Wr  at  height  of  burst  hR,  one  would  expect  at  a  range  of  Rr  a  reflected 
peak  overpressure  in  the  Mach  reflection  region  equal  to  that  obtained  from  a  free-air  burst  with 
a  yield  of  RpWj. 

An  indication  of  the  reflection  factor  to  be  expected  in  the  Mach  reflection  was  obtained 
from  the  information  presented  in  Figure  HI.  2  relating  the  reflection  factor  for  TNT  to  the  scaled 
height  of  buret  for  a  particular  (and  not  necessarily  ideal)  reflecting  surface.  This  and  the  additional 
information  for  pentolite  and  the  related  data  for  nuclear  bursts  indicated  that  scaled  heights  of 
bursts  could  be  selected  such  that  the  reflection  factor  in  the  far  Mach  region  would  be  between  0,5 
and  2.  0,  approximately.  In  Figure  III.  1  are  shown  the  R  W  curves  for  these  two  extremes  relating 
the  reflected  peak  overpressure  with  the  scaled  slant  range.  The  curves  in  Figure  III.  3  provide  the 
same  information  for  a  nonscaled  range. 
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FIGURE  UI.  1.  PEAK  OVERPRESSURE  VS  SCALED  DISTANCE  FOR 
SPHERICAL  TNT  IN  FREE-AIR  AT  SEA  LEVEL 
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mm, *L£.cr/o*/  sac  to  a  -  mf , 
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FIGURE  III.  2.  REFLECTION  FACTORS  FOR  TNT  EXPLOSIONS 
AT  LOW  HOB  OVER  VARIOUS  SURFACES 
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FIGURE  ni.  3.  PEAK  OVERPRESSURE  VERSUS  RANGE  FOR  TNT 
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Since  the  above  method  of  predicting  the  reflected  peak  overpressure  was  limited  to  the  far 
Mach  region  and  in  order  to  avoid  having  the  missile  model  experience  two  pressure  peaks  in  the 
regular  reflection  region,  the  charge  sise,  height  of  burst  and  the  ground  range  were  selected  such 
that  the  triple  point  would  be  above  the  uppermost  portion  of  any  model,  regardless  of  its  orientation. 

The  information  in  the  literature  provided  the  necessary  relation  between  the  Mach  stem 
height  and  the  range  for  various  heights  of  burst.  From  prior  qualitative  tests,  it  was  known  that  the 
missile  models  would  sustain  damage  from  a  30-lb  charge  of  TNT  at  a  ground  range  of  15  to  26  feet, 
depending  on  the  orientation  and  construction  of  the  model  itself.  The  model  positioned  with  its  lon¬ 
gitudinal  axis  perpendicular  to  ground  represented  the  most  severe,  triple-point  height  requirement. 
For  this  arrangement  and  a  ground  range  of  12  feet,  the  scaled  height  of  hurst  would  be  limited  to  a 
scaled  HOB  of  0.  8  or  less  or,  for  a  charge  weight  of  30  lb  of  TNT,  a  burst  elevation  above  ground  of 
2-1/2  feet  or  less.  Referring  to  Figure  III.  2,  this  height  of  burst  would  provide  a  reflection  factor 
of  approximately  1. 8  and  therefore  ground  reflection,  peak  overpressures  at  the  model  of  from 
approximately  15  to  100  psi  depending  on  the  ground  range. 

It  was  recognized  that  such  factors  as  the  actual  ground  reflection  characteristics  at  the  test 
site,  the  type  of  charge  and  the  charge  shape  would  all  somewhat  influence  the  blast  phenomena. 
Accordingly,  pressure  measuring  devices  were  used  to  determine,  among  other  things,  the  extent 
of  the  deviation  (if  any)  of  the  predicted  peak  overpressure.  (See  Appendix  II.  ) 

The  results  of  these  pressure  readings  (shown  as  points  superimposed  on  the  curves  shown 
in  Figures  III.  1  and IH.  2)  are  especially  encouraging  particularly  in  view  of  the  uncertainties 
normally  associated  with  blast  phenomena  measurements.  It  is  seen  that  with  the  exception  of  the 
30- lb,  H4  shot  and  the  5-lb,  El  and  E2  shots,  a  ground  reflection  factor  of  between  1 . 6  and  1 . 7  for 
the  30- lb  shots  and  of  0.  7  for  the  5-lb  shots  is  quite  valid.  The  values  for  Rf  were  used  for  predict¬ 
ing  the  peak  reflected  overpressures  in  those  tests  where  pressure  measurements  were  not  obtained. 
Table  Ill.  1  summarizes  the  pertinent  peak  overpressure  and  reflection  factor  data. 

In  addition  to  the  overpressures,  the  other  blast  parameters  which  significantly  influence  the 
extent  to  which  the  model  is  damaged  are  the  positive  phase  duration  and  the  positive  impulse.  The 
time  of  arrival  of  the  blast  wave  is  of  importance  in  the  establishment  and  interpretation  of  the 
model's  high-speed  photography  data. 

.Available  free-air  curves  for  arrival  time,  positive  phase  duration  and  positive  impulse  for 
TNT  provided  suitable  approximations  for  adjusting  the  pressure  recording  and  high-speed  photo¬ 
graphy  equipment.  Table  III.  2  offers  a  comparison  between  measured  and  predicted  data.  A  sum¬ 
mary  of  the  pressure  data  is  given  in  Table  III.  3. 

III.  2.  HIGH-SPEED  PHOTOGRAPHY  DATA 

The  variation  of  the  radial  displacement  at  x  =  L/2,  <f>  =  0*  with  time  (where  t  =  0  denotes  the 
arrival  of  the  blast  wave  at  the  cylinder's  leading  edge)  is  shown  in  Figures  III.  4  through  III.  13. 

Each  of  the  curves  was  fitted  to  points  obtained  from  either  prints  or  projections  of  each  individual 
frame . 


The  curves  in  Figures  III.  8  and  III.  10  need  some  additional  interpretation.  The  model  (CB2) 
used  in  shot  HI  had  been  previously  subjected  to  four,  5-lb  explosions  in  shot  series  G.  These 
smaller  (but  repeated  exposures)  resulted  in  a  slight  permanent  set  at  the  cylinder's  top,  leading 
edge  adjacent  to  the  grooved  (for  simple  support)  end  plates.  In  the  subsequent  shot  HI,  this  edge 
of  the  shell  moved  out  of  the  groove  at  approximately  1 . 7  msec  (see  Figure  III.  8)  after  the  arrival 
of  the  blast  wave.  In  so  doing,  the  subsequent  response  of  the  shell  was  no  longer  that  of  a  cylinder 
with  simple  supports  at  each  boundary.  However,  prior  to  that  time,  the  shell's  response  followed 
the  theoretically  predicted  displacement  pattern. 

The  same  mode  of  failure  was  experienced  by  Model  CB3  (Figure  III.  10)  in  shot  H4.  Here, 
the  shell's  leading  edge  at  the  top  end-plate  began  to  move  out  of  the  groove  at  approximately  1. 5 
msec.  Although  tue  shell  did  subsequently  return  to  zero  displacement  at  x  =  L/2,  $  =  0*,  the  condi¬ 
tion  of  simple  support  was  no  longer  in  effect.  For  this  reason,  the  simple- support  response  was 
limited  to  that  portion  of  the  displacement  curve  prior  to  1. 5  msec. 
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TABLE  III.  1.  PEAK  OVERPRESSURE  IN  MACH  REFLECTION 
REGION  (EXPERIMENTAL  DATA) 


Ground^  ^ 


Shot 

No. 

Charge^ 

(lb) 

HOB<2> 

(ft) 

Range 

(ft) 

(A) 

RF  =  0.  5 

RF  =  2-0 

(B) 

_R£_ 

(C) 

(D) 

(E) 

El 

5 

0.  33 

10 

12.  8 

33.  0 

1.  32 

25.  1 

>5.0 

2.0 

E2 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

29.  3 

1.  6 

E3 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

18.8 

0.9 

FI 

30 

1.  5 

20 

10.  5 

27.  2 

1.  89 

27.  0 

23.0 

1.  6 

F  2 

30 

1.  5 

16 

16.  3 

47.0 

1.  89 

46.4 

41.0 

1.7 

G1 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

19.  7 

0.9 

G2 

5 

0.  33 

10 

12.  8 

33.  0 

1.  32 

25.  1 

16.  1 

0.  7 

G3 

5 

0.  33 

10 

12.  8 

33.  0 

1.  32 

25.  1 

16.  5 

0.7 

G4 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

N.  A. 

N.  A. 

H2 

30 

2.  5 

20 

10.  5 

27.2 

1.  89 

27.  0 

24.  6 

0.  7 

H3 

30 

2.  5 

20 

10.  5 

27.  2 

1. 89 

27.0 

23.0 

1.  6 

H4 

30 

2.  5 

15 

18.  5 

55.9 

1.  89 

54.0 

53.0 

1.6 

J1 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

13.  1 

0.  6 

J2 

5 

0.  33 

10 

12.  8 

33.0 

1.  32 

25.  1 

16.2 

0.  7 

J3 

5 

0.  33 

5 

39.  2 

140.  0 

1.  32 

102.0 

64.0 

0.  7 

N.  A.  Not  available. 

(1)  Flaked  TNT. 

(2)  Height  of  burst. 

(3)  Distance  from  ground  zero  to  model  and/or  pressure  pickup. 

(A)  Estimated  range  of  peak  overpressure  in  Mach  reflection  region  (psi); 
Kirkwood -Brinkley  theory  -  see  Figures  III.  1  and  III.  3. 

(B)  Reflection  factor  for  hard  clay  surface  as  function  of  HOB  (see  Figure  III.  2). 

(C)  Peak  overpressure  based  on  Rp  from  Column  (B). 

(D)  Experimentally  measured  peak  overpressure  (psi). 

(E)  Adjusted  reflection  factor  based  on  measured  peak  overpressure. 
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TABLE  HI.  2.  TIME  OF  ARRIVAL,  POSITIVE  PHASE  DURATION 
IN  MACH  REFLECTION  REGION  (EXPERIMENTAL  DATA) 


1 

f 

Shot 

Charge*1* 

hob*2* 

Ground*2* 

(A) 

(B) 

(C) 

(D) 

1 

1 

* 

No. 

(lb) 

Range  (ft) 

(msec) 

(msec) 

(msec) 

(msec) 

D1 

30 

2.5 

26 

11.5 

N.  A. 

_ 

_ 

t 

D2 

30 

2.5 

22 

9.3 

7.9 

— 

— 

9 

> 

( 

D3 

30 

2.5 

16 

5.0 

5.3 

— 

— 

l 

El 

5 

0.25 

10 

3.4 

N.  A. 

2.6 

2.7 

t 

E2 

5 

0.25 

10 

3.4 

N.  A. 

2.6 

2.8 

\ 

\ 

E3 

5 

0.  25 

10 

3.4 

N.  A. 

2.6 

2.2 

f 

FI 

30 

1.5 

20 

7.4 

N.  A. 

5. 1 

4.2 

* 

F2 

30 

1.5 

16 

4.9 

N.  A, 

4.0 

4.3 

f 

G1 

5 

0.25 

10 

3.4 

4.5 

2.6 

2.4 

} 

G2 

5 

0.  25 

10 

3.4 

3.9 

2.6 

2.3 

G3 

5 

0.25 

10 

3.4 

3.9 

2.6 

2.5 

G4 

5 

0.25 

10 

3.4 

3.8 

2.6 

2.7 

i 

HI 

30 

2.5 

20 

7.4 

7. 1 

5. 1 

N.  A. 

H2 

30 

2.5 

20 

7.4 

N.  A. 

5. 1 

6.0 

H3 

30 

2.5 

20 

7.4 

N.  A. 

5. 1 

W4.  5 

» 

H4 

30 

2.5 

15 

3.7 

N.  A. 

3.9 

N.  A. 

i 

i 

J1 

5 

0.25 

10 

3.4 

N.  A. 

2.6 

3. 1 

J2 

5 

0.25 

10 

3.4 

N.  A. 

2.6 

3.4 

J 

J3 

5 

0.25 

5 

1.2 

N.  A. 

1.6 

1.1 

N.  A.  Not  available. 

(1)  Flaked  TNT. 

(2)  Height  of  burst. 

(3)  Distance  from  ground  zero  to  model  and/or  pressure  pickup. 

(A)  Estimated  time  of  arrival  -  ta 

(B)  Recorded  time  of  arrival  -  ta  (from  high  speed  photography). 

(C)  Estimated  positive  phase  duration  -  t+ 

(D)  Recorded  positive  phase  duration  -  t+ 


TABLE  III.  3.  SUMMARY  OF  PEAK  INCIDENT  OVERPRESSURES 
AND  POSITIVE  PHASE  DURATIONS 


(1) 

E:  end-on  loading;  S;  side-on 

(4) 

E  -  elastic;  F  -  fixed; 

loading;  O;  oblique  loading 

S  -  simple  support 

(2) 

Peak  incident  overpressure 

(5) 

Experimental  -  See 
Table  III.  1  or  III.  2 

0) 

Positive  phase  duration 

Shot 

Model  (1) 

pi  (2) 

t+  (3) 

Boundary 

Model 

No. 

Orientation 

(psi) 

(msec) 

Restraint  (4) 

CC1 

A1 

£ 

4 

7 

E 

A2 

1 

10 

6 

A3 

1 

f 

24 

4 

B1 

S 

18 

5 

B2 

23 

4.5 

r 

B3 

32 

4 

CC2 

A1 

2 

9 

A2 

6 

7 

A3 

17 

5 

B1 

29 

4 

B2 

38 

3 

CC3 

Cl 

31 

5 

CC4 

Cl 

£ 

7 

22 

6 

CC4 

C2 

£ 

27 

5 

CC5 

Cl 

O 

13 

7 

CC5 

C2 

O 

17 

6 

CA1 

Cl 

E 

22 

6 

F 

CA1 

C2 

E 

27 

4 

CAZ 

Cl 

s 

> 

27 

4 

CA3 

D1 

15 

6.5 

D2 

19 

5.5 

D3 

42 

3.8 

CB1 

FI 

27 

4 

S 

CB1 

F2 

46 

4 

CB2 

G1 

20  (5) 

2.4(5) 

G2 

16  (=) 

2.3  (5) 

G3 

17  (5) 

2.3  (5) 

G4 

18 

2.5 

HI 

25 

5. 1 

CB3 

H3 

23  (5) 

4.5 

CB3 

H4 

53  (5) 

3.9 

PI 

J1 

- 

25 

3 

P2 

J2 

- 

25 

3 

1 

P2 

J3 

- 
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FIGURE  III.  4  RADIAL  DISPLACEMENTS  AT  x  =  L/.L  4.  =  0'  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  III  S.  RADIAL  DISPLACEMENTS  AT  x  =  L /2,  4>  =  0’  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  III  6  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  <»  -  0°  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  HI.  7.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  +  =  0*  OBTAINED  FROM  HIGH 
SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  in.  8.  RADIAL  DISPLACEMENTS  AT  X  *  L/2,  ♦  <=  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  III.  9.  RADIAL  DISPLACEMENTS  AT  x  ■  L/2,  ♦  *  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  UI.  10.  RADIAL  DISPLACEMENTS  AT  x  *  L/2,  *  *  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 


V 

'L 


* 

I 

U 

o! 

VI 

«** 

Q 


5 

Q 


FIGURE  III.  11.  RADIAL  DISPLACEMENTS  AT  x  *  L/2,  $  =  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  IU.  12.  RADIAL  DISPLACEMENTS  AT  x  =  L/2,  4>  =  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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FIGURE  UI.  13  RADIAL  DISPLACEMENTS  AT  x  =  L/2.  *  =  0*  OBTAINED  FROM 
HIGH  SPEED  PHOTOGRAPH  OF  CYLINDRICAL  MODELS 
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